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ABSTRACT 


An experimental investigation of a large-radius 
pinch discharge has revealed that the cylindrical cur- 
rent sheet develops a broad and diffuse anode attachment 
which, under certain conditions, divides into two dis- 
tinct regions of high current density. Furthermore, the 
existence of precursor current conduction at downstream 
anode locations was inferred from data which showed a 
substantial increase in arc voltage in response to insu- 
lation of the central portions of the pinch anode. Other 
experiments indicated that both of the above noted anode 
phenomena occur in a parallel-plate accelerator, and also 
that the analogous cathode processes arise at somewhat 
later times in the discharge history. Experimental study 
of various anode processes in a quasi-steady MPD arc has 
clarified certain aspects of the current conduction and 
power loss mechanisms in these devices. Measurements of 
anode fall -voltage and local anode current density at 
"matched" operating conditions have defined an overall 
inverse relationship between these two quantities which 
is reciprocal in form. This result translates into a sig- 
nificant decrease in fractional anode power loss for quasi- 
steady MPD arc operation in the megawatt range of pulse 
power. For matched and overfed mass flows, anode current 
conduction in this device is accomplished by random thermal 
flux of electrons from the adjacent arc plasma; for under- 
fed conditions, substantial electric fields develop around 
the anode in order to support anode current conduction at 
this low plasma density situation. 
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OiAPTER I 
INTRODUCTION 

By detailed mission analysis, and by consideration of 
the equation of motion for a rocket, it has been shown that _ 
interplanetary space missions can only be effectively carried 
out by space thrusters which provide specific impulses of over 
1000 sec; i.e., propellant exhaust velocities in excess of 
10,000 m/sec. Exotic chemical and nuclear rockets have been 
extended to specific impulses in the 800 sec range, but it is 
becoming increasingly apparent that a different propulsion 
concept is required to achieve those values of specific im- 
pulse demanded by future space missions. The acceleration of 
gases by electrical heating and/or by electric and magnetic 
body forces has demonstrated a very high specific impulse ca- 
pability at relatively low thrust densities, and now defines 
a promising new propulsion system known as electric space pro- 
pulsion. The penalty of low thrust capability is not severe 
in the weak gravitational fields of interplanetary space; and 
is more than counterbalanced by the much higher ratio of pay- 
load mass to initial rocket mass achievable by means of the 
high exhaust velocities of electric space thrusters. 

Jahn^ subdivides the field of electric propulsion 
into three classes: 

1. Electrothermal propulsion, wherein the pro- 
pellant gas is heated electrically, then expanded in a suitable 
nozzle. 

2. Electrostatic propulsion, wherein the propel- 
lant is accelerated by direct application of electric body 
forces to ionized particles. 
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3. Electromagnetic propulsion, vdierein an ionized 
propellant stream is accelerated by interactions of external 
and internal magnetic fields with electric currents driven 
through the stream. 

Both electromagnetic as well ,as certain electrothermal 
plasma accelerators consist basically of a high-current arc 
discharge passing through a streaming plasma between two suit- 
ably designed and insulated electrodes. The thermal arc jet, 
operating on the electrothermal principle, takes advantage of 
very high arc column temperatures to increase the enthalpy of 
the propellant gas, and then converts this enthalpy to stream- 
ing energy by ejjpanding the gas through a nozzle which also 
functions as the anode. At low mass feed rates and high arc 
currents, the propellant flow becomes almost totally ionized 
and the magnetic fields become sufficiently strong to enable 
the electromagnetic mode of plasma acceleration to take over, 
causing significant increase in specific impulse. Devices with 
coaxial electrode configuration, in which the interaction of 
high arc currents with self- or imposed-magiietic fields directly 
accelerates . an ionized flow of plasma, became known as magneto- 
plasmadynamic (MPD) arc jets. 

In general the electromagnetic class of plasma accel- 
erators can be further subdivided into those that operate in 
the steady or continuous mode and those that operate in the 
unsteady or pulsed mode. Steady accelerators, of which the MPD 
arc jet has demonstrated the best performance, employ a dc power 
supply capable of delivering at least 3000 A of continuous cur- 
rent at kilowatt power levels; but, are restricted to operation 
in this low power range because of severe heat transfer to the 
electrode and limitations of the' dc power supply. Unsteady 
plasma accelerators normally have pulse times on the order of 

4 

p^secs with current levels above 10 A and employ a variety of 
electrode configurations? parallel rail, parallel plate, co- 
axial, and cylindrical z-pinch. The acceleration process be- 
gins when the voltage of a capacitor network is suddenly trans- 
ferred to the accelerator electrodes causing the ambient gas- 
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fill to ionize and conduct current in a narrow region at the 
position of minimum inductance. The narrow .current -distri- 
bution, or current sheet, the shape of which depends on the 
accelerator geometry, is then acted upon by the j x B mag- 
netic interaction forces which drive the sheet along the elec- 
trode channel, forcing it to entrain and accelerate the ambient 
gas to interesting propulsion velocities. Because of the very 
short duration of the driving-current pulse, devices of this 
type do not suffer from serious electrode heating and hence 
can operate at extremely high-power levels which appear to 
entail high specific impulses and high thermal efficiencies. 
However, pulsed accelerators are limited to a thrust efficiency 

of about 50% because of the inelastic piston behavior of the 
2 

current sheet. Furthermore, the extremely short pulse times 
place severe demands on the gas-feeding and power -switching 
systems if the accelerator is to be effectively operated in 
a repetitive manner. 

3 

Recent studies have demonstrated the feasibility and 
advantages of utilizing the so-called quasi-steady mode of 
accelerator • operation, wherein self-field MPD arcjets are 
operated in a pulsed manner but essentially steady, high-power 
operation is achieved over 100 psec to several msec pulse du- 
rations. The overall concept involves repetitive pulsing of 
the accelerator at repetition rates commensurate with available 
space power supplies, thereby affording the opportunity for 
variable thrust accomplished via simple duty cycle adjustment. 

These devices, known as quasi-steady MPD arcjets, promise high 

' 4 

thermal efficiencies and high specific impulses, both inherent 

to the high power regime; as well as high thrust densities along 
with tolerable heat transfer and- mean power demands. Much of 
the work reported here was carried out on a quasi-steady MPD 
arc which is described briefly in Chap. Ill and with more de- 
tail in references 3 and 5. 

Considerable experimental work has been done on ther- 
mal and MPD arcjets to assess the thermal efficiencies of these 
devices. . The results indicate that thermal efficiencies 
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range from 10 to 70%, with most of the power loss appearing 
at the anode, and also that these efficiencies tend to improve 
at high power levels. The investigations also indicate that 
most of the anode heat transfer is due to the flux of current- 
carrying electrons which deposit energy gained in the anode 
fall sheath at the anode. Hence, the magnitude of the anode 
fall voltages which arise in steady and quasi-steady MPD arc 
discharges determine, to a large extent, the amount of anode 
power loss, and thus the thermal efficiency of these devices. 
Part of this thesis, therefore, is devoted to an investigation 
of the anode voltage drop in a guasi-steady MPD arc. 

Much of the' early work on the electrode regions of the 
arc discharge^^ concentrated on the cathode, as it was con- 
sidered to be the electrode most critical for the maintenance 
of the discharge. Furthermore,' most of the arc discharge 
studies reported in the literature have been carried out on 
so-called classical arcs, distinguishable by currents- in the 
ampere range and pressures near one atmosphere, conditions far 
removed from those of interest to electric propulsion. In re- 
cent years, studies of the anode region have been extended to 
thermal arcjets and steady, low-:power MPD arcs where the intro- 
duction of streaming plasmas, strong magnetic fields, and co- 
axial geometry has added considerable complexity to the prob- 
12—14 

lem. Extension of MPD arc .operation into the multimega- 

watt power range, through the application of the quasi-steady 
concept, incorporates a host of new physical processes at the 
anode region which cannot be inferred from studies on low- 
power MPD arcs. Thus, in order to further the understanding 
of the microscopic and macroscopic processes prevailing at the 
anode of quasi-steady, high-power MPD arcs, in particular those 
of current conduction and heat transfer, this experimental re- 
search program was undertaken on the quasi- steady MPD device 
at Princeton. 

Interest in anode phenomena in this laboratory grew 
originally from the early experiments of Black on an 8-in. 

5 

pinch discharge employing currents of over 10 A for jusec pulse 
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durations. Among other contributions, these studies pointed 
out the occurrence of a peculiar current sheet bifurcation; 
i.e., an apparent splitting of the current sheet at the mid- 
way point of its rapid collapse towards the chamber center. 

Later experiments revealed that the anomalous behavior was, 
in fact, related to a broad, diffuse, and sometimes divided 
current sheet anode attachment . Further investigations of . 
this phenomena uncovered a precursor conduction current at 
the anode surface which bears a strong influence on the re- 
sistive voltage of the pinch discharge. These studies were 
later extended to a pulsed parallel-plate accelerator • where 
several modifications in electrode phenomena were observed. 

The results of these investigations are presented in Chap. IV 
and V of this thesis. 

The experiments on the pinch and parallel-plate de- 
vices were concluded when it became evident that better under- 
standing of the current sheet anode attachment and anode pre- 
cursor conduction would require detailed investigations of 
various plasma instabilities and surface phenomena, somewhat 
beyond the scope of this work. At the same time, it was de- 
cided that the previously described quasi-steady MPD arc would 
be a more interesting facility in which to study anode pro- 
cesses. The electrode geometry, along with the quasi-steady 
nature of this MPD discharge, facilitate' detailed study of 
the anode region with electrostatic and magnetic probes; while 
the present interest in quasi-steady plasma acceleration fur- 
ther emphasizes the importance of extending anode studies to 
this device. As a result, the focus of this thesis shifts after 
Chap. V to experiments on and analysis of the anode region of 
the quasi-steady MPD arc. Thus,, the present investigation of 
anode phenomena includes experimental study at various dis- 
charge geometries: cylindrical pinch, linear parallel-plate, 

coaxial; and at two different operation conditions: unsteady 

and quasi-steady. 
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CHAPTER II 

REVIEW OF ANODE PROCESSES IN ARC DISCHARGES 

A. INTRODUCTION 

The term^ electric arc, is used tp define a particular 
kind of electrical discharge characterized by high current and 
very low impedance. Electric arc discharges have been studied 
extensively for many years, and have recently become relevant 
to the concept of plasma acceleration because of their ability 
to conduct very high currents at reasonably low voltages, while 
producing very high plasma temperatures. High arc currents 
along with their attendant strong magnetic fields can be used 
to accelerate ionized gases to very high velocities by means 
of the magnetic interaction, or j x B body forces. Another 
class of accelerators, the electrothermal variety, capitalize 
on the high temperatures produced in the arc column, up to 
50,000 °K, to increase the enthalpy of the propellant gas which 
is subsequently' expanded through a nozzle. This chapter will 
acquaint the reader with some of the physical processes which 
participate at the anode-plasma junction of conventional elec- 
tric arcs; and will also examine several more recent studies 
carried out on the thermal arcjets and MPD arcs referred to in 
Chap . I . 

B. DESCRIPTION OP THE CLASSIC ELECTRIC ARC 

The electric arc is distinguished from other types of 
gas discharges such as the corona, abnormal glow, and normal 
glow, by the fact that the cathode is sufficiently heated by 
ion bombardment to enable the thermionic emission of electrons 
to occur, and this mechanism permits the arc discharge to attain 
its very low impedance. 
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The classical electric arc can be roughly broken down 
into three regions, each of which plays an important and dis- 
tinct role essential to the overall physical process: 

1) A positive column , which extends over most 

of the arc length and whose purpose is the conduction of cur- 
rent through the main body of the gas. The electric fields 
are rather weak in this part of the arc, and 'the plasma, usually 
in thermal equilibrium, consists of a strongly radiating mix- 
ture of electrons, ions, and neutral atoms at nearly the same 
temperature, from 5000 to 50,000 having corresponding 

degrees of ionization from a few percent to essentially 100 
percent. The random thermal velocities of the electrons in 
this region far exceed their mean drift velocities in the weak 
field along the column, and thus the predominant ionization 
mechanisms are thermal electron-atom collisions and photoion- 
ization, Current conduction is mainly via a diffusion-domi- 
nated drift of the electrons. 

2) A cathode fall , which exists at the cathode 

surface and whose function is the transfer of current across 

the gas-cathode junction. This requires high electric fields 

relative to those in the column and results in a positive 

space-charge accumulation. The cathode fall is typically less 

than 20 V and the cathode surface emits electrons prolifically 
3 7 2 

(10 to 10 A/cm ) by a combination of- thermionic, photoelec- 
tric, and field emission processes. Only in this part of the 
arc do ions contribute significantly to the current conduction 
during their acceleration through the cathode fall. 

3) An anode fall , which exists at the anode sur- 
face and whose function is the transfer of current across the 
gas-anode junction, also requiring high electric fields. The 
anode fall is roughly of the same magnitude and spatial extent 
as the cathode fall (on the order of the dominant collision 
mean free path) , and the principle ionization processes involve 
the field-accelerated electrons. In general, the anode does 
not emit ions, and current conduction is almost entirely by 
electrons which leave the positive column via a diffusion- 
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dominated drift and are accelerated towards the anode by the 
anode fall. On the other hand, ions move from the anode to- 
wards the cathode due to the electric field and hence a neg- 
ative space charge develops. To maintain this steady flow 
of ions required for electrical neutrality of the column, and 
to guarantee that the spatial extent of the space-charge re- 
gion be small, ions, must be produced in a thin sheath contig- . 

19 

uous to the anode surface. According to Hocker and Bez the 
ion production in this sheath may occur either by field ioni- 
zation, wherein the electrons fall freely through the anode 
fall space and ionize neutral atoms at collision only after 
they have acquired an energy at least as great as the ioniza- 
tion energy of the gas or vapor, or by thermal ionization, 
wherein the electrons undergo several collisions in the anode 
fall so that their directed energy becomes thermalized and 
only the electrons in the high energy "tail" of the energy 
distribution engage in ionizing collisions with neutrals. The 
value of the anode fall voltage is then determined by the 
energy requirements for ion production. 

A major difference between the cathode fall and the 
anode fall lies in the amount of charge carrier production re- 
quired in the two regions. It can be shown that Ij^/® ions and 
electrons per second must be generated at the anode fall whereas 
I^/e electrons per second must be produced at the cathode, where 

I . and I are the ion and electron components of the total arc 
1 e 17 

current in the quasi-neutral plasma of the column. Moreover, 

near the anode the current is almost entirely electronic and 
hence only about I. /I (proportional to the ratio of the rela- 
tive drift velocities in the column of the two species, about 
10 ionizations per incident electron are needed. Hence, 
almost’ all of the current-carrying electrons reach the anode 
without undergoing inelastic collisions and deliver the energy 
they have gained in the fall directly to the anode. This 
results in a large arc power loss at the anode (see Sec. E) , 
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C. PIASMA COOLING AT THE ANODE 

X8 

According to Bez and Hdcker/ the thickness of the 
anode fall sheath is of the order of the electron mean free 
path in the anode-adjacent plasma. This calculated thick- 
ness agrees with experimental measurements of Block and 
23 

Pinkelnburg who found, from moving probe measurements, a 

sheath thickness of 2 p for the carbon arc at atmospheric 

pressure. Adjacent to this fall sheath, a second region 

exists in which the field intensity is lower than in the fall 

but still a factor of approximately 50 to 100 higher than in 

23 

the column. Experimental evidence shows that this region 
extends at most to a distance € =0.1 mm from the anode sur- 
face. It is caused by the cooling effect of the anode sur- 
face which decreases the electrical conductivity in this zone; 
it will be referred to as the C - layer. 

) The temperature dependence of the electrical conduc- 

tivity of argon at a pressure of 1- atm,' as calculated by Busz 

24 

and Finkelnburg, is shown in Fig. 2-1. Two regions exist in 
the CT - T variation, one in which <5" rises very rapidly with 
temperature, and a second in which CT" is roughly constant for 
increasing plasma temperature. An approximate value of the 
temperature separating these two regions is T 10, 000 °K, 
which is also very close to the temperature of the plasma ait 
the outer edge of the - layer ( y = 6 ) . Obviously at 

y > ^ the temperature drop due to the cooling of the plasma 
by the anode has no effect upon the electrical conductivity, and 
hence the field intensity in the outer region of the thermal 
boundary layer ( S'-r > y > ^ ) is comparable to that in 

the arc column; i.e., relatively low. In the region 0<y<^^ 
however, where the electrical conductivity decreases rapidly 
with decreasing wall distance and thus decreasing plasma tem- 
perature, the electric field intensity must increase sharply 
in order to maintain current continuity. Thermal equilibrium 
between electrons and heavy gas particles no longer exists at 
this high field intensity because. of the strong energetic 
coupling between the electric field and the electrons, and the 
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TEMPERATURE x 10’ °K 


TEMPERATURE DEPENDENCE OF ELECTRICAL 
CONDUCTIVITY AT A PRESSURE OF I otm. 
{AFTER BUSZ-PEUCKERT a FINKELNBURG (24)) 


FIGURE 2- I 
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weak coupling between the heavy particles and both the field 

and the electrons. Thus/ the electrons receive additional 

energy from the field which they cannot transfer to the heavy 

particles, resulting in an increase in the electron temperature 

above that of the heavy particles T^, in the € “• layer. A 

■qualitative representation of the distribution of temperature, 

25 

voltage, and electric field is given in Fig.- 2-2 / where the 
thermal boundary layer thickness, S'x , is used as a ref- 
erence distance. The anode fall thickness is too small to be 
shown on this scale. Henceforth the combination of the anode 
fall and the £ -layer shall be referred to as the anode volt- 
age drop. Electrostatic probe measurements in the pinch dis- 
charge and quasi-steady MPD arc, presented in Chap. IV and VI, 
were taken at distances from the anode which are too great to 
allow for a distinction between the anode fall and the £ -layer. 
Hence, the voltage measured at the anode surface will include 
contributions from both of these regions, and therefore can only 
be used to define the anode voltage drop rather than the anode 
fall voltage specifically. 

D. ANODE PALL THEORIES 

The only detailed theoretical investigation of the anode 

fall region is that developed in a series of papers by-Bez and 
^0 22 

Hocker. They cited three main functions of the anode 

fall: 

1) The production of ions which then diffuse 
through the cross section of the column towards the cathode. 

2) The directed velocity of the ions, accelerated 
towards the column by the anode fall, must be converted into a 
thermal velocity distribution corresponding to the thermal tem- 
perature of the positive column or to that of the adjacent 

£ -layer. 

3) The accomplishment of the transition from the 
relatively low temperature anode to the high temperature of the 
positive column or adjacent £ -layer; to this end the fall 
probably heats the cold gas adjacent to the anode surface. 
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TEMPERATURES , VOLTAGE AND ELECTRIC 
FIELD INTENSITY NEAR THE ANODE SURFACE 

(REF.25) 


FIGURE 2-2 
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It is Tcnown that^ with equal energy^ the differ- 
ential coefficient of ionization of electrons is appreciably 
larger than that of the ions. For this reason/ only elec- 
trons are taken into consideration in the ionization processes. 
Bez and Hocker have distinguished two different types of ioni- 
zation by electrons in the fall region; field ionization and 
thermal ionization. Their characteristics .are discussed below. 


1. Field Ionization 

The term, field ionization, applies to the state of 
affairs in which electrons fall freely through the anode fall 
space using the energy thus acquired to ionize gas or vapor 
atoms at collision. To investigate this possibility, Bez and 
Hocker divided the. anode space-charge zone into two different 
model zones, an "acceleration zone" and an "ionization zone." 
The "acceleration zone" is that part of the anode fall zone in 
which the voltage, starting from the column, grows to the ioni- 
zation voltage of the gas, V^, while the "ionization zone" is 
described by the region in which the voltage grows from the 
ionization voltage to the anode fall voltage, V^. For fhe 
special case of a carbon arc burning freely in air at atmos- 
pheric pressure, the authors calculate the extent of these 

—3 —3 

areas as 0.4 x 10 cm and 1.8 x 10 cm, respectively. 

As Bez and Hocker consider the extent of the anode 
fall area, for this type of ionization, as less than one elec- 
tron mean free path, they calculate the motion of the electrons 
in the acceleration zone in accordance with the laws of free 
fall. If one further ignores the thermal energy at entry into 
the acceleration zone, then the electrons will have at point x 
the energy corresponding to the potential V(x) . 

To describe the ionization zone, from X. = X© to 
X ~ Xa t the authors use the following expression for the 
ion production over a distance dx 




X" 


(2-1) 
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where Cj C '/r ) is an approximation of the' coeffi- 

cient of differential ionization, and is the electron 

mean free path for ionization. This expression can be in- 
tegrated to give an exponential form, and since cIx/'Xq 
the exponential can be approximated by a linear form 
yielding 



It is assumed that. the ion current vanishes at the anode sur- 
face, X=Xa ' above expression then gives 



By ignoring the ion current density and the counter 

18 

field diffusion of ions toward the anode, Bez and Hdcker 
approximate the voltage V(x) in both the acceleration and 
ionization zones by a "generalized Langmuir formula" 


VM = ^ 2 ,. 3 


3 rm<s. TT ^ -/e 


'/3 


2 e 



(2-4) 


in which Q, is an adjustable parameter used to fit the ex- 
perimental anode fall data. Application of the Langmuir for- 
mula, derived from a collisionless model, to the ionization 
zone is based on the assumption that this region is much 
smaller than an electron collision mean, free path, and thus 
most electrons fall. freely through the, ionization zone to 
the anode surface. 

Equation (2-3) can now be readily integrated by using 
the above expression for the voltage variation in the fall 
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zone. The result is 



3 /MaTT^ ^ Xe 

2& J Cx 
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where the authors use Cx^o) ^ -^e and ^><o) 

which are respectively the electron and ion current densities 

a * • 

at the anode end of the positive column. When /e and jt 
are written in terms of the electron and ion mobilities, the 
following expression for \/a , as an implicit function of T, 
is derived 


( 



7T 3^'^e^'^qf'^E \ ^ 
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A e A 4 a. 
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where T is the temperature of the plasma at the anode' end of 
the column (assumed to be in local thermodynamic equilibrium) , 
and Xto. and are the ion-atom and electron-heavy par- 

ticle collision mean free paths. 

By adjusting the value of (X, , Eg. (2-6) can be 

used to predict the correct anode fall voltage, = 20 V, 
for the low current carbon arc burning in nitrogen at atmos- 
pheric pressure with a gas temperature of about 6000 °K. How- 
ever, the authors noted that the calculated anode fall volt- 
ages rise sharply at temperatures beyond about 9000 °K. This 
results mainly because of the rapid . decrease of / due 

to increasing electron-ion collisions, as the degree of ioni- 
zation begins to rise (see Fig. 2-3). In fact, Eq. (2-6) pre- 
dicts an anode fall voltage of 35 V for the above described 
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arc when the .nitrogen temperature is raised to 10,000 
Since this result contradicts the experimental data, Bez and 
HocTcer claim, therefore, that the field ionization model can- 
not describe the ion production in the anode fall zone for 
arc column temperature beyond about '9,000 °K. Instead, ion 
production at the anode fall of high temperature arcs proceeds 
via the thermal ionization mechanism which 'is described next. 

Since Joule-heating of the column plasma depends upon 
the square of the current density ( , it follows that 
high arc temperatures are normally associated with high-cur- 
rent density arcs. Therefore, the distinction between field 
and thermal ionization in the anode fall can also be made on 
the basis of current density. For example, the low current 
carbon arc, with an anode current density of 40 A/cm and 
anode fall of 20 V, exhibits the properties of field ioniza- 
tion at the anode, whereas the high-current carbon arc, with 

an anode fall 8 to 10 V lower, is characterized by the thermal 

17 

lonxzation mechanism. 

2. Thermal Ionization 

The term, thermal ionization, denotes the state of 
affairs in which electrons undergo several collisions in the 
anode fall so that their velocities are randomized and ap- 
proximately thermal with ah electron temperature, J'eA / 
somewhat above that of the heavy particles. There are suffi- 
cient electrons in the high energy "tail" of the energy dis- 
tribution to produce the necessary ionizing collisions and 
resulting positive ions; hence, the anode fall need not be 
as large as when field ionization prevails, where it must be 
at least as large as the first ionization potential of the 
vapor or gas. The anode fall thickness will be greater than 
for field ionization because the electrons must now make 
several collisions in this region, and therefore the fall 
zone must be larger than 

A simple energy balance for the thermal production 
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of ions takes the form 

^*Va - ( S’/z MTe/\) (2-7) 

where /v'a is the power density deposited in the anode fall, 
^ I is the ion current density leaving the cathode end of the 
fall, and the last term represents the power delivered to the 
anode surface by the current-carrying electrons, which have 
become thermalized in the fall at a temperature ~T~ef\ * Prom 
Eg. (2-7), the following estimate of the anode fall voltage 
requirement can be derived, 

\4 = \ca Vt -h 5- ykTsA (2-8> 

^ (rmtJ Xet ^ ~e~ 

where and have again been expressed in terms of 

the relative mobilities of the electrons and ions in the 
column plasma. For this type of ionization, the mean thermal 
electron energy, JiTefi , can be much smaller than the 
ionization energy of the gas or vapor; in addition, the ratio 
(cff^e / first term small even at high gas 

temperatures where the ratio Xia/Xet becomes large (see 
Pig. 2-3) . Thus, Eg. (2-8) demonstrates that the anode fall 
voltage requirement for ion production by the thermal ioniza- 
tion process can, indeed, be considerably less than the ioni- 
zation potential of the gas. 

E. ANODE ENERGY BALANCE 

An anode heat transfer model for steady plasma ac- 

. 9 

celerators, formulated by Shih, Pfender, Ibele, and Eckert, 
includes a variety of transfer processes such as energy trans- 
fer by electrons carrying the current, heat transfer by com- 
bined conduction and convection, and by radiation. The 
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electrons carry to the anode surface an amount of energy^ 

5 kT^/2 per electron, corresponding to their thermal energy 
in the neighboring column plasma, plus the amount of energy 
gained through the anode fall, V^. In addition, an amount 
of energy proportional to the work function of the anode ma- 
terial, , appears also as heat input. There has recently 

been some discussion about whether the contribution due to 
the electron thermal energy in the column should be 5 bT^/2 
or 2 This question is taken up in the Appendix with 

the conclusion that 5 kT^/2 is the better representation of 
this effect. The convective heat transfer from the plasma 

to the anode surface can be expressed in terms of the enthalpy 

27 

difference between the plasma stream and the anode surface. 

The total heat input to the anode is then given by 


- he h/\( is - -t Pra -t (p^) ^2-9) 


where 


ht - 
is. = 

u = 

Aa = 

PRA = 
X = 

Te = 


heat transfer coefficient at the anode 

surface based on enthalpies 

enthalpy of the gas outside the boundary 

r 

layer 

enthalpy which the gas assumes in the 
immediate proximity of the anode surface 
surface area of the anode 
radiative heat transfer from the cathode 
and neighboring plasma to the anode 
total arc current 

electron temperature at the anode end of 
the column 


The total arc current is used in this expression 
rather than the electron current specifically because, at the 
anode, the arc current is carried almost entirely by electrons. 
Furthermore, as was noted in Sec. B, only about 0.1 to 1.0 
percent of electrons entering the anode fall give up their 
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energy in ionizing collisions with neutral atoms, therefore 
this energy need not be subtracted from the electron contri- 
bution to the anode energy flux. 

A confirmation of the proposed anode heat transfer 

model was obtained experimentally by measuring the total heat 

9 12 

transfer to the anode for various arc currents, ' with the 
result that plots of Pa vs J' gave a straight Tine, in- 
dicative of constant values of ~Ye Va • was con- 

cluded from this data that the dominant contribution to the 
anode energy loss is due to the electron current, particularly 
at high arc currents, and hence the equation for anode heat 
transfer can be approximated by 

P f\ ^ ^ -''kjk + \/ f\ + (P/i ^ (2- ID) 

The expression for heat transfer to the cathode con- 

' T * 9 

tains a similar term which is linear in the ion current, . 
Since for high-current arcs J~ >> CJj / and since electron 
emission at the cathode is an endothermic process, it is evi- 
dent that the cathode power loss is much less significant than 
anode power loss. There is also considerable empirical evi- 

9 

dence supporting this claim. Hence, it follows that most of 
the attempts to improve the thermal efficiency of plasma ac- 
celerators have been focused on alleviating the anode power 
loss . 

Much experimental work has been devoted to the mea- 
surement of anode heat losses in high-current arc jets and MPD 
arcs. It has been reported that the relative heat losses at 
the anode depend on arc parameters such as current, mass flow 
rate, and applied magnetic field. Without detailing the de- 
pendence of the anode losses on these quantities, the average 

results quoted by a number of investigators will now be given. 

6 

Schoeck found that the total heat flow to the anode for a 
high pressure thermal arcjet, operating between 50 and 150 A 
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at a pressure of 1 atm, was between 75 and 90% of the total 

7 

arc power ♦ Schneiderman and Patrick, as well as Ducat i, 

8 

Muehlberger, and Treat, also working with thermal arcjets, 
observed that as much as 60 to 70% of the input power was 
lost to the anode. Shih et al. reported anode power losses 
from 30 to 60% in argon and amonnia-fed MPD arcs. Finally, 
Nerheim and Kelly, in a review of MPD arc technology, found 
that the more technically promising engines had a median range 
of thermal efficiency of around 50 to 70%, with most of the 
power losses appearing at the anode - 

p. Transpiration cooling of the anode 

In view of the substantial heat load which must be 
carried by the anode of high -current plasma accelerators, 
attempts were made both to cool. the anode and to lower the 
net anode heat loss by a process known as transpiration 
cooling of the anode; i.e. , the injection or forcing of a 
suitable coolant through an anode which has been fabricated 
from a porous material. The transpiring gas readily pro- 
vides the cooling needed to prevent ablation of the anode and, 
in addition, feeds a large portion of the energy .normally lost 

to the coolant back into the arc region where it may serve a 

25 

useful purpose. Schoeck et al. injected argon gas through 
a porous graphite anode back into their high intensity argon 
arc, and observed that, although between 50 and 80% of the arc 
input power was transferred to the anode, up to 80% was re- 
coverable by this process. In fact, thermal efficiencies as 
high as 80% were attainable. Pfender et al. employed a con- 
stricted, wall-stabilized argon arc, with power inputs from 
2.5 to 16 kW, to compare a water-cooled copper anode and a 
porous graphite transpiration-cooled anode. They found that 
the anode heat losses were reduced by a factor of two, from 
20% to only 10% of the total input power, when transpiration 
cooling was used. Hence, it would appear -that transpiration 
cooling of the anode is an effective means of improving the 
thermal efficiency of low power arc thrusters. However, as 
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will be shown in Chap. VI, high-power, quasi-steady MPD arcs 
promise high thermal efficiencies, in the 80% range, without 
the complications of this anode cooling procedure. 

G. ANODE PRECURSOR PHENOMENA IN PULSED ACCELERATORS 

29 

Lie, Ali, and Chang have observed precursor effects 
along the anode surface ahead of the main current sheet in 
both parallel-plate and coaxial accelerators- They used mono- 
chrometers to detect preionization in the ambient gas ahead 
of the current sheet, and observed that the level of preioni- 
zation was much higher at the anode than at the cathode for 
helium, nitrogen, and argon prefilling pressures of 0.05 to 
0.5 Torr. In addition, they observed that this preionized 
gas near the anode carried a considerable fraction of the dis- 
charge current, particularly when the filling pressures were 
low. It was concluded that the anode preionization and pre- 
cursor conduction are related neither to the initial high- 
current discharge near the back insulator, nor to the radiation 
emitted by the main discharge current, but rather, are a 
vestige of the breakdown conduction processes which occur 
along the entire electrode surface. 

The existence of this anode surface current was used 
by Lie et al. to explain the incline of the current sheet to- 
wards the anode in parallel-plate and rail accelerators. This 
so-called, current sheet tilt, wherein the anode attachment 

runs ahead of the main portion of the current sheet, has been 

30—3 2 

well documented in previous works, and was also observed 

in the present experiments (see Chap. V) . Since the current 
sheet near the anode has a forward axial current component, 
due to precursor surface conduction, this part of the sheet 
experiences a j x B force directed toward the cathode, and 
this tendency amplifies itself as the current sheet travels 
down the electrode gap resulting in a tilt of the sheet as 
a whole. 

33 

It is instructive to note that Lovberg has observed 
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a planar/ untilted current sheet in the case of a hydrogen- 
filled (0.3 Torr) parallel-plate accelerator; and also that 
Lie et al. found no appreciable preionization and no axial 
current along the anode of their accelerator when hydrogen 
was used. These two observations suggest that the phenomena 
of current sheet tilt and precursor anode conduction are in- 
deed causally related as proposed by Lie et al. . In addi- 
tion/ it can be postulated that these two related phenomena 
are somehow associated with the large ion Larmor-radii of the 
heavier gases, such as nitrogen and argon, which in some cases 
approach inter-electrode dimensions. Hydrogen, on the other 
hand, which produces planar current sheets and no appreciable 
pre ionization, indicates ion Larmor-radii at least an order 
of magnitude smaller. 

Other experiments, performed by Lovberg on a parallel- 
plate accelerator with hydrogen-filling, showed that a 100 V 

anode sheath was present, with the cathode sheath being much 

33 

smaller. This large anode sheath voltage was explained in 

terms of a current sheet model, in which ions are presumed 

to be the major current carriers, and for which Lovberg found 

32 

supporting ejqjerimental evidence. The small Larmor-radii 
predicted for the hydrogen ions support this conclusion, since 
the hydrogen ions will be therefore able to make complete rev- 
olutions within the current sheet, and thus may drift towards 
the cathode, carrying current in the process. With ions 
carrying the current and moving away from the anode, a large 
excess of electrons will be left there, and the resulting neg- 
ative space charge causes the large voltage drop at the anode 
surface. The rather small cathode sheath voltage arises be- 
cause the current sheet is an effective ionizing agent, thus 
making electron emission at the cathode and ionization in the 
cathode fall unnecessary. 

34 

Kislov et al. have mapped the potential contours in 
a quasi-steady coaxial plasma accelerator, operating at cur- 
rent levels from 20 to 50 kA and nitrogen mass flow rates from 
1 to 3 g/sec. They also observed large potential drops at the 
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anode surface (30 to 110 V) , or about 40% or more of the 

total discharge voltage, while the drop near the cathode was 

relatively small. Working with a similar device, KDvrov et 
35 

al. measured the magnetic fields and current dxstrxbutions, 

and discovered that the discharge current pattern appeared to 

"slide" along the anode surface and thus extend well down- 

34 

stream of the cathode current zone. Kislov et al. attempted 
to explain the sizable anode voltage drop by suggesting that 
the electrons are prevented from reaching the anode by virtue 
of their large Hall parameters, which enable them to undergo 
a cross-field drift parallel to the anode surface in response 
to the electric and magnetic fields, which are respectively 
perpendicular and parallel to the anode. Hence, rather than 
being absorbed at the anode, the electrons form a negative 
space-charge layer there, resulting in a substantial anode 
voltage drop. Furthermore, it was proposed by Kovrov et al. 
that the electron Hall drift parallel to the anode in the 
downstream direction causes the current distribution to ex- 
tend along the forward portion of the anode, the so-called 
"anodic sliding." The measured values of electron tempera- 
ture and density at the anode of this device confirm that the 
electron Hall parameters are indeed large enough for this 
effect to occur. 

34 

The conclusions of Kislov et al., which relate the 

anode voltage drop in their accelerator to the large electron 

Hall parameters, correlate well with earlier experiments of 

3 6 

Karkosok and Hoffman, in which electrode voltage drops in 
a MHD channel, with a current of only 1 A and a pressure of 
1 atm, were measured. Their results indicate that the anode 
voltage drop doubled, from 20 V to 40 V, for an increase in 
the electron Hall parameter from very small values to about 
2.5. 

37 

More recently, Dammann et al. have observed a pre- 
cursor in a thin layer near to the anode in front of the con- 
vergent cylindrical shock wave of z-pinch discharges in 
helium and argon. The precursors, which were investigated 
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by interferometric methods, appear to separate from the 
shock wave at higher filling pressures, and then run ahead 
along the anode toward the chamber center. The authors have 
attempted to explain their observations in terms of gasdy- 
namic effects. 

H. REMARKS 

In this chapter a discussion of plasma cooling at 
the anode and a theory of the anode fall zone have been pre- 
sented- The physical understanding developed in these sec- 
tions will be employed in later chapters to explain certain 
experimental results pertaining to the anode region of pinch 
and quasi-steady MPD discharges. In addition, the anode heat 
transfer model is applied to the quasi-steady MPD arc in 
order to compute the anode power loss and thermal efficiency. 
Finally, the discussion of the anode precursor phenomena aids 
our attempts to understand and explain similar phenomena in 
our -pulsed accelerators. 
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CHAPTER III 

.EXPERIMENTAL FACILITIES AND DIAGNOSTIC TECHNIQUES 

A. INTRODUCTION 

The study of anode phenomena in high-current dis- 
charges was carried out, for reasons stated previously, under 
a variety of discharge geometries and operating conditions. 

For example, pulsed or unsteady discharges, employing a prop- 
agating current sheet as a plasma acceleration mechanism, were 
studied in both cylindrical z-pinch and linear, parallel-plate 
geometries at current levels from 100 to 300 kA and pulse times 
from 3.5 to 20 psec. The quasi-steady mode of electromagnetic 
plasma acceleration was investigated in a coaxial device em- 
ploying current levels from 4.4 to 42 kA and pulse durations 
from 60 to 600 psec. Argon was used as the operating gas for 
each of these experiments because it is relatively inexpensive 
and because its relatively high molecular weight affords the 
opportunity for high thrust density. 

B. CYLINDRICAL Z-PINCH FACILITY 

The large-radius linear pinch discharge is charac- . 
terized by the early formation of a uniform cylindrical current 
sheet at the outer chamber radii, the configuration of lowest 
inductance, followed by the rapid collapse of this current 
structure towards the chamber center as it is driven by self- 

«jk 28 

field, j X B forces. ' The main components of the discharge 

apparatus are the cylindrical discharge chamber, a gas- triggered 

39 . ... 

switch, and a series of high-voltage capacitors distributed 

15 

along a transmission line. Briefly, the discharge chamber is 
8-in. in diameter and 2-in. high; the electrodes, of 3/4-in. 
aluminum, are separated by a 1/4-in. -thick, 2-in; high Pyrex 
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sidewall (Fig. 3-1) . The switch unit consists of a similar 
8-in. pinch chamber fitted with a gas-feeding Plexiglas in- 
sert. The discharge is initiated by injecting argon into the 
previously evacuated switch chamber, thereby causing the switch 
to break down and transfer the capacitor bank potential to the 
electrodes of the prefilled main discharge chamber. The cham- 
ber electrodes are both held at ground potential during the 
charging of the capacitor bank by means of a 1 kn. ballast 
resistor. The electrical energy is initially stored in a bank 
of 20 capacitor units, each with nominal 2.5-pf capacitance 

and 15-nH inductance, rated to 10 kV. The capacitors are ar- 

. 15 

ranged along an LC transmxssion line. The driving-current 
waveform can be altered by simply changing the distribution 
of the capacitor units along the line. The discharge chamber 
is normally prefilled to 100 p argon; and the entire system is 
evacuated to several p. after each discharge by means of a me- 
chanical pump and then flushed with argon gas. Chamber pres- 
sures are monitored with a standard Pirani gauge. The dis- 
charge current is measured with a Rogowski coil, of standard 
40 . . 

design connected to a passive RC integrator circuit. Data 
were generally recorded photographically from a Tektronix Type 
555, dual-beam oscilloscope. 

C. PARALLEL-PLATE ACCELERATOR 

The parallel-plate experiments are, in many respects, 

similar to the cylindrical pinch; the only major differences 

are in the geometry and dimensions of the apparatus (Pig. 3-2) , 

In the present arrangement 3/4-in. aluminum electrodes, 6-in. 

wide and 48-in. long, are separated and supported by 2-in. high 
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Plexiglas sidewalls. A rectangular gas-triggered switch is 
employed to initiate the discharge sequence. Electrical energy 
is stored in a capacitor bank consisting of 40 3.2-pP capac- 

itors rated at 30 kV. The four rows of 10 capacitors each can 
be arranged in either parallel or series to vary the magnitude 
and shape of the driving-current waveform; from 120 kA x 20 jusec 
to 4.4 kA x 600 ^sec. The discharge is again initiated when 
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the switch is injected with argon, causing it to break down 
and transfer the 10 kV bank potential to the accelerator elec- 
trodes. A uniform, planar current sheet is initially estab- 
lished at the back insulating wall, the position of lowest 
inductance; and then propagates down the length of the elec- 
trode channel, driven by self-magnetic forces. The total dis- 
charge current and accelerator width are such that the re- 
sulting current densities are sufficiently high to guarantee, 

at least in the early phases of the discharge, a well-defined 
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current- sheet. In some ej^perxments, nylon insulation was in- 
laid along the electrodes in order to arrest the propagating 
current sheet and induce a stabilized current pattern. 

The supporting apparata employed with this facility 
are similar to those associated with the 8-in. pinch device; 
a mechanical pump to evacuate the system after each discharge, 
a Pirani gauge, a Rogowski coil, and a Type 555 dual-beam os- 
cilloscope. The accelerator channel is normally prefilled to 
100 /a argon. This device affords several diagnostic advantages 
over the 8-in. pinch facility; longer streamwise dimensions, 
longer time scales, and easier photographic observation. The 
disadvantages are; a less reproducible discharge and greater 
three-dimensionality of the current sheet, both caused by the 
side-wall distortion of the magnetic field lines. 

D. QUASI -STEADY MPD ARC FACILITY 

The quasi-steady mode of arc operation affords distinct 
diagnostic advantages over the unsteady operation of the 8-in. 
pinch and parallel-plate facilities- The pulse duration is 
sufficiently long that the troublesome transient and inductive 
effects can subside, yet is sufficiently short to permit the 
application of simple unsteady diagnostic techniques. In 
addition, problems of probe destruction, encountered in the 
hostile environments of steady-flow MPD arcs, are not serious 
in'' the quasi-steady operating regime. The other attractive 
features of the quasi-steady MPD arc have already been dis- 
cussed in Chap. I, and therefore will not be restated here. 
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The quasi-steady MPD arc facility, employed in the 
present investigation of anode phenomena, is described in de- 
tail in reference 5. This device closely replicates the con- 
ventional magnetoplasmadynaraic arc configuration, but on a 
scale one' order of magnitude larger to permit detailed in- 
terior diagnostics. A schematic diagram of this apparatus is 
shown in Pig. 3-3. Briefly, it consists of a cylindrical dis- 
charge chamber with a 3/4-in. -diameter conical tungsten cath- 
ode, and an aluminum anode with a 4-in. -diameter orifice, driven 
by a 130 fiF x 10 kV capacitor line, arranged to provide a vari- 
ety of rectangular current pulses ranging from 4.4 kA x 600 psec 
to 140 kA X 20 psec. Propellant gas, normally argon, is in- 
jected through six calibrated orifices in the arc chamber end- 
plate from a high pressure reservoir, abruptly established by 
the end wall stagnation of the flow in a simple shock tube. 

An adjustable bleed line from the same shock tube triggers the 
discharge circuit switch, similar to the one described in 
Chap. Ill, Sec, B, thereby permitting correlation of the mass 
flow pattern with the arc current profile. The arc chamber 

and thruster assembly are mounted in a 3-ft. diameter x 6-ft. 
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long Plexiglas vacuum tank, which is evacuated down to pres- 

_5 

sures of about 10 Torr by a 6-in. -diameter oil diffusion pump. 

The operation sequence of the system begins with the 
pressurization of the shock tube driver section to 45 psig, 
normally, and the charging of the capacitor line to 10 kV. 

The shock tube diaphragm is ruptured, and the pressure wave 
is first bled into the evacuated switch and then establishes 
a pressure reservoir at the end wall of the driven section, 
which in turn feeds the argon gas into the chamber through the 
six end-wall orifices. The bleed-line length is adjusted so 
that the switch breaks down and transfers the 10 -kV bank po- 
tential to the accelerator electrodes, only after steady mass 
flow has been established in the arc chamber. The steady 
phase of operation, when both current and mass flow patterns 
have become stabilized, can be achieved over a large portion 
of the driving-current pulse; all arc data are collected during 
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this steady phase. The experiments on this quasi-steady MPD 
arc facility are carried out at currents from 4.4 to 42.0 kA 
and argon mass flow rates from 1.2 to 36.0 g/sec. 

E. MAGNETIC PROBE DIAGNOSTICS 

40 

Magnet xc probes, or magnetic induction coils, have 
been employed to measure the magnetic field patterns in each 
of the previously described experimental facilities. These 
probes consist of small, multi-turn coils of insulated wire, 
which generate a voltage across their leads in response to a 
change in the magnetic flux threading the coil- .The probe 
signal, by Faraday’s law, is proportional to the r ate of change 
of magnetic field at the sensor location. Integration of the 
probe signal by a passive RC integrator circuit provides a 
record of the magnetic field over the discharge history. By 
assuming azimuthal symmetry of the discharge. Ampere’s law can 
be used to relate the local magnetic field strength to the total 
current enclosed within a radius defined by the coil location 
in the discharge chamber. Hence, 
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(3-1) 


where S is the plane circular surface defined by rotating the 
radius, of the probe location through 360°. Thus 

B -c/i = (3-2) 

I ^ 

where Tp is the radius o’f the probe location; and assuming 

i 

azimuthal symmetry, this integral yields 
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which relates the local magnetic field strength, B, at a ra- 
dius r , to the total current closed by that radius, J , 

p ^ enclosed 

Thus, contours of equal enclosed current can readily be plotted 
from a knowledge of the magnetic field distribution; these con- 
tours can then be used to deduce the distribution of current 
density in the arc discharge. 

Several turns of thin Formvar wire, wound on- a 1/8-in.- 
diameter coil, provide sufficiently large signals and suffi- 
ciently fine spatial resolution to accurately measure the strong 
magnetic fields of the high-current, pinch and parallel-plate 
devices. However, in order to measure the weak magnetic fields 
at the anode of the quasi-steady MPD arc, much more sensitive 
probes are required. The probes used for this purpose consist 
of 100 turns of #40 Formvar wire wound on a 3/16-in. -diameter 
nylon spool.- The probes are bent to allow access to the back 
portions of the anode; and the signals are fed to a Type 555 
dual-beam oscilloscope through an integrator circuit with 
RC = 1500 jasec, a time one order of magnitude longer than the- 
current ■ pulse in order to avoid excessive droop in ' the inte- 
grated signal. 

F. INNER VOLTAGE DIVIDER 

Voltage studies were carried out on both- the cylindrical 
pinch and parallel-plate facilities, wherein the resistive volt- 
age drop across the current sheet must be separated from the 
inductive voltage contribution associated with the rapidly prop- 
agating current distribution. The pure resistive component of 
the discharge voltage can be monitored by the proper placement 
of an inner divider voltage tap; i.e., an insulated metal rod 
which completes a circuit between the two electrodes such that 
no magnetic flux is enclosed. In the pinch geometry the volt- 
age top must be located at the chamber center, whereas a volt- 
age top inserted at the downstream end of the parallel-plate ac- 
celerator serves the same purpose. Magnetic fields fold out be- 
hind the propagating current sheet; i.e., interior to the cur- 
rent sheet and the return conductor path. Hence, it can readily 



35 


be seen that a measuring circuit consisting of the current 
sheet, the two electrodes, and a downstream voltage hap, will 
not enclose any magnetic flux (see Pig. 3-4) . A Tektronix 
Model 6013 voltage probe is used to monitor the total voltage 
drop of this measuring circuit, which is simply the resistive 
voltage across the current' sheet neglecting the voltage drops 
along the electrodes and metal voltage tap. 

G. TWIN LANGMUIR PROBE TECHNIQUE 

Twin Langmuir probes have been successfully employed 
to measure the electron temperature and number density at sev- 
eral .locations inside the arc chamber, and adjacent to the 
anode of the quasi-steady MPD facility. These probes consist 
of two separately shielded, and biased, 1/8-in. long exposed 
tips of 3-mil tungsten wire, oriented parallel to the flow axis 
to minimize streaming effects, and placed sufficiently close 
together to sense essentially the same plasma. The probe ra- 
dius is less than the anticipated electron-ion mean free path 
in the region to be probed and larger than the Debye length, 

in order to minimize the effects of particle collisions and 
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sheath size on the measurements. 

The twin probe technique involves the attachment of 
different resistor pairs to the two probe tips, and the mon- 
itoring of the voltage developed across these resistors with 
Tektronix P6006 voltage probes (Pig. 3-5) . A differential 
Type G pre-amplifier is used with a Tektronix Type 555 dual- 
beam oscilloscope to record the voltage difference between the 
two probe tips. In this way, a single discharge of the accel- 
erator will provide sufficient information for the determination 
of local electron temperature, provided that both probes oper- 
ate in the electron retarding portion of the Langmuir probe 
characteristic. This simple calculation is based on the fol- 
lowing reasoning. The electron current to a probe tip, which 
is biased to retard electron flux, is given by 

Iq = .e A p ^nrie Ce^ ^ e Vp ^ 


(3-4) 
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SCHEMATIC OF TWIN LANGMUIR PROBE 
AND CIRCUITRY 


FIGURE 3-5 
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where A is the area of the probe tip, n and T are the 
p ^ e e _ 

electron number density and temperature respectively, C is 
the mean thermal electron speed, and is the probe poten- 
tial relative to plasma potential. When both probe tips 
operate on the electron retarding branch of the characteristic, 
the ratio of the relative electron currents is 



where \/f and V 2 refer to the actual probe tip voltages since 
the plasma potential value cancels out of the difference. The 
advantage of the twin electrostatic probe technique rests with 
the fact that l 2 » Vj^, and are measured at the same time 
in the same discharge, thus assuring that n and C will can- 
cel precisely in the above ratio, and thereby eliminating the 
problem of shot-to-shot fluctuations in electron density. Fur- 
thermore, it is required that the probe tips draw a sufficiently 
large electron current that the ion current to the probe can 
be neglected, therefore 


I 


e f 




cy 




(3-6) 


where 2 actual currents drawn by the respective 

prpbe tips, calculated from the measured probe voltages, 

..and known resistor values, 2 * Thus, taking the logarithm 
of Eg. (3-5) and using the above approximations yields 



e 
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or 

JlTe = ® ev < 3 - 7 ) 

To reiterate/ this expression for electron temperature 
is valid only when both probes are retarding the electron flux 
and the electron flux to the probe still greatly exceeds the 
ion flvix. Probe data must be collected for a number of dif- 
ferent resistor pairs in order to establish when the proper 
domain of probe operation has been achieved, and when the 
probes cease to retard electrons and pass over onto the elec- 
tron-collecting branch of the probe characteristic. This tran- 
sition occurs when the probe potential becomes equal to the 
plasma potential and results in a noticeable change in probe 
behavior. At this special condition, the probe current is 
given by 

eAp(m^CeAf-) 



The electron number density, n^, is thus readily obtained from 

the previously determined electron temperature, T^, the probe 

surface area, A , and the measured probe current at the tran- 

sition from the electron-retarding to the electron-collection 
* 

regions, I . 

Figure 3-6 displays a typical plot of twin Langmuir 
probe data for 17.5 kA x 5.9 g/sec operation of the quasi-steady 
MPD arc; the twin probe was located at a chamber radius of 
3/4 in., and axial position 1 in. upstream of the front anode 
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face. The quantity Q ^ i /Xa ^ is plotted 

against the larger of the two probe currents for several dif- 
ferent resistor pairs. Several features of this plot are im- 
mediately evident. The plateau between 1.75 and 2.75 A corre- 
sponds to probe operation in the high-current portion of the 
electron retarding branch. In this region of the character- 
istic, therefore the electron 

temperature is read directly from the ordinate to be 

kT = 1.35 + 0.15 eV 
e — 

The rise in the curve at around 3.0 A indicates that 

the probe has reached the transition region and hence has a 

potential equal to that of the surrounding plasma. At this 

condition the electron density is given by Eq. (3-8) . Sub- 
* 

stituting I =3.0+0.5A from the above plot, and inserting 
T^ = 1.35 + 0.15 eV into Eq. (3-8), yields the following value 

n = 1.3 + 0.4 X 10^® 
e — 

The departure from the plateau at the low currents 
reflects probe operation on the low current portion of the 
electron- retarding branch where the approximation, Ig I. 
is violated. The quantity no 

longer equals the local electron temperature and begins to 
decrease as the two probe potentials start to converge. 
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CHAPTER IV 

EXPERIMENTS ON THE 8-INCH PINCH DEWCE 
A. BACKGROUND 

Early experimental work carried out on the 8-in. pinch 
device with associated pulse-forming network, described in 
Chap. Ill, yielded certain anomalous data on current sheet dy- 
namics which eventually led to the present studies of anode 

15 

phenomena. These early studies compared the observed cur- 
rent sheet trajectories with those computed from the so-called 
"snowplow model" of the pinch event. This model considers the 
current sheet to be an impermeable front, driven by self-mag- 
netic fields, which entrains and imparts substantial velocity 
to the gas as it sweeps towards the chamber center. The term 
"inelastic piston" aptly describes this acceleration model. 

For most situations considered, the agreement between theory 
and experiment was rather good, except in certain cases where 
the current sheet appeared to divide spontaneously at about 
mid-chamber radius into two concentric sheets of like sign, 
which then followed each other to the center of the chamber. 
This peculiar event was first observed when a relatively flat 
current pulse of 200 kA and 5.5 psec duration was used to drive 
an 8-in. pinch discharge into 100 p of ambient argon fill. So- 
called "current sheet bifurcation" was detected by means of 
unintegrated magnetic probe records which departed from the 
normal single peak configuration, representing the passage of 
a single current sheet, and began to display two distinct 
maxima, implying the passage of two similar current sheets. 
After the inception of current sheet bifurcation, the slower 
of the two signal peaks continued to follow the predicted 
snowplow trajectory, whereas the leading peak ran foirward at 
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approximately twice the snowplow velocity. Figure 4-1 dis- 
plays the trajectory plots for the bifur eating- sheet along 
with the theoretical snowplow trajectory; also shown is the 
driving-current waveform used in this experiment. It was 
noted that a factor of approximately two separates the leading 
sheet velocity from the snowplow velocity. This suggests an 
"elastic piston" model, wherein particles are elastically re- 
flected off of the moving current sheet, thereby acquiring a 
velocity which is twice that of the sheet. 

B. MAGNETIC PROBE STUDIES 

In order to test the elastic piston hypothesis, fur- 
ther magnetic probing of these interesting pinch discharges 
was carried out. The above quoted data were taken with mag- 
netic probes constrained to the midplane of the pinch device; 
to explore the axial development of the bifurcation event, 
magnetic probe readings were taken at several different axial 
positions in the chamber. Figure 4-2 displays the magnetic 
probe responses at a chamber radius of 1 in. along with the 
400 kA X 3.6 psec driving- current pulse, also found to produce 
anomalous sheet behavior. With the probe located 1/4 in. from 
the anode. Fig. 4-2a, two distinct maxima are seen, implying 
the passage over the sensing coil of two distinct zones of 
high-current density. In Pigs. 4- 2b and 4-2c, the probe is 
moved from a position 1/2 in. from the anode to one at the 
chamber midplane, with the result that the bS/''d't (unin- 
tegrated probe signal) maxima have moved closer together and 
appear to coalesce into a single peak similar to that observed 
at a probe position 1/2 in. from the cathode (Fig. 4-2d) . This 
single maxima class of response is maintained everywhere from 
this location up to the cathode surface. 

Figure 4-2 yields two other important pieces of in- 
formation. The front-running peak near the anode 

is ahead of the midplane response, thus agreeing with the pre- 
viously observed tilt of the current sheet wherein the anode 
attachment runs ahead of the main portion of the sheet, see 
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Chap. II, Sec. G. Also, the magnitude of the maximum probe 
response near the anode is about one— half that of the probe 
near the cathode, and more broadly spread, implying a more 
diffuse current density structure at the anode. 

A collection of data of this type now suggests a sub- 
stantially different image of the bifurcation process. Namely, 
that the current extends to the cathode in a single, well-de- 
fined current sheet, but enters the anode in a more diffuse, 
double-humped attachment. The two current density maxima are 
widely separated at the anode surface, and coalesce to the 
single sheet structure at some "triple point" away from the 
surface thus forming a "Y" -shaped configuration. Fig. 4-3. 

As this structure propagates toward the center of the chamber, 
the separation of the anode attachment widens and the "triple 
point" progresses axially away from the anode, steadily eroding 
more and more of the original single sheet. Magnetic probes 
placed on the midplane thus would observe a single sheet at 
large radii, and would be subjected to the well-developed anode 
bifurcation at small radii, recall Sec. A. 

Further studies with magnetic probes have demonstrated 
that the phenomena of anode bifurcation can occur for a va- 
riety of driving-current waveforms and ambient argon chamber 
pressures. The event is favored by relatively slow rise time 
pulses of either flat or double-humped shape, and was observed 
over a broad range of chamber pressure. However, the extent 
of the anode attachment structure and the rate at which it 
erodes the current sheet depend on both of these discharge 
parameters. Furthermore, the anode bifurcation event is pe- 
culiar to large-radius pinch discharges with relatively long 
discharge times. Burton, working with a 5-in. -diameter 
pinch discharge driven by a 300 kA x 1. 2 /nsec current pulse, 
observed only well-defined, reproducible current sheets with 
no evidence of a broad and divided anode structure. On the 
other hand, Kolpin^^ observed double maxima in current density 
in a large-radius inverse pinch, employing 16-in. -diameter 
electrodes, a 600 kA x 14 jisec current pulse, and ambient 
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helium at 500 p. From these observations it would appear 
that anode bifurcation is also a function of discharge time 
and geometry. The 1.2 /nsec pinch time of Burton's small 
radius discharge is simply too short for this anode structure 
to develop to any observable degree. 

C. PHOTOGRAPHIC STUDIES OF ANODE BIFURCATION 

Following the initial observations of this unusual 
current sheet behavior, attempts were made to map the current 
density distribution in the "anode foot" region, but this was 
precluded by an inordinate scatter in the data, a significant 
result in itself for these otherwise highly reproducible 
closed-chamber discharges. The unusual irreproducibility of 
these measurements is interpreted to indicate that the dis- 
charge current experiences considerable difficulty in estab- 
lishing a viable conduction mechanism near and into the anode 
surface, and thus is vulnerable to unspecified instabilities 
and turbulence in this region. 

In deference to the above-mentioned irregularity of 
the process, a systematic photographic study was undertaken 
to catalog the visual aspects of the domain before attempting 
further accumulation of electrical data. For this purpose a 
series of radial-view, Kerr-cell photographs were taken of the 
discharge. This was accomplished by removing a section of 
the outer return conductor, allowing radial optical access to 
the interior of the chamber through the glass wall insulation 
(Fig. 4-4) . All of the displayed photographs were taken with 
a 100 p argon discharge in a 8-in. diameter chamber driven by 
a 5.5 psec x 200 kA current pulse. A magnetic probe, visible 
in the photographs, was also inserted axially into the cham- 
ber to provide d8/^t records along with Kerr-cell 

photographs. 

The first series of photographs, displayed in Fig. 4-5, 
shows the progress of the pinching current sheet from 1.4 psec 
to 3.8 psec after initiation of the discharge. These pictures 
again show that the luminous front does divide and diffuse 
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near the anode, and that the disturbance grows towards the 
cathode as the current sheet propagates inward. 

Correlation of such Kerr-cell photographs with cor- 
responding magnetic probe records allows one to relate the 
luminosity front to the current conduction regions of the 
sheet. Figure 4-6 compares photographs with the appropriate 
magnetic probe records for three different positions of the 
probe. In Fig. 4-6a, both the luminosity pattern and the 
magnetic probe data indicate a single current conduction zone 
at the midplane. For Fig. 4-6b the magnetic probe has been 
positioned 1/2 in. from the anode. The photograph shows the 
probe standing between the two "legs" of the bifurcated cur- 
rent sheet. The corresponding probe trace shows that one 
current conduction region has swept over the probe prior to 
1.8 psec and that a second current zone is about to pass the 
probe. Apparently, then, each part of the bifurcated lumi- 
nosity front has a current conduction region associated with 
it. Figure 4-6c, with the magnetic probe at a position 1/4 in. 
from the anode, further supports this conclusion. 

In the course of this qualitative photographic survey, 
it seems reasonable to try to establish the pressure depend- 
ence of the anode foot processes. Figure 4-7 shows photographs 
of current sheets propagating into ambient argon at pressures 
of 50 p, 100 p, and 200 p at approximately the same radial 
positions. In each case the driving current is constant at 
about 200 kA for the times indicated. Each photograph shows 
the familiar current sheet bifurcation near the anode, with 
little drastic change in scale over the range tested. In 
fact, the degree of similarity among the three figures would 
seem to indicate that chamber pressure is not a critical fac- 
tor in determining the current sheet behavior near the anode. 

Under certain conditions another aspect of the anode 
process has been observed photographically. When a highly 
polished aluminum electrode is employed, a distinct fila- 
mentary "spoking" is found to occur at a random spot on the 
anode. Figure 4-8 illustrates the situation: at the early 
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times the familiar bifurcated anode attachment is observed, 
but later, at about 2.6 fiisec, a luminous "spot" is seen to 
form on the anode surface. This luminous spot intensifies 
as the current sheet propagates inward. This is not to be 
confused with the central chamber luminosity, seen in the 
early photographs, which is caused by a process known as 
"bleed through, " in which some of the intense luminosity pro- 
duced when the sheet collapses at the chamber center is able 
to penetrate the nitrobenzene shutter of the Kerr cell and 
register on the film. From many photographs of this type, it 
is concluded that the "spoke" does not form at any one par- 
ticular point on the anode, but seems to attach at a different 
point for each discharge. Little is presently known about the 
current densities associated with this "anode spoke" but it 
seems reasonable, however, that such instabilities give rise 
to the irreproducible magnetic probe data mentioned earlier. 
Finally, this "spoking" phenomena may be interpreted to mean 
that the discharge current is having difficulty establishing 
a feasible conduction mechanism at the anode, and must there- 
fore resort to unusual and erratic means of making the tran- 
sition across the metal-plasma interface. 

D. ARC VOLTAGE MEASUREMENTS 

Further experiments associated with the study of cur- 
rent conduction properties at the anode of the 8-in. pinch 
discharge have focused on measurements of the arc voltage in 
this device. Since unsteady discharges of this kind have 
very large inductive voltage contributions, it is necessary 
to employ an inner divider voltage probe (see Chap. Ill, 

Sec. F) to measure the pure resistive component of the total 
discharge voltage, referred to as the arc voltage. Since the 
arc voltage is a good index of the energy delivered to the 
electrodes and to Joule-heating of the arc plasma, and since 
the broad and diffuse anode attachment of the pinching cur- 
rent sheet probably involves certain dissipation of energy, 
it seems reasonable to seek, now, some correlation between the 
arc voltage and the current distribution near the anode. 
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For this purpose, the 8-in. pinch discharge was mon- 
itored with an inner voltage divider constructed of an insu- 
lated metal rod inserted across the chamber electrodes to 
complete a measuring circuit of which the current sheet is 
part. The signals were presented to the oscilloscope via a 

g 

Tektronix Model 6013 voltage probe of 10 ohm impedance. The 
relevance of the current sheet anode attachment, as well as 
other features of anode current conduction, in determining 
the resistive voltage drop, or arc voltage, was assessed by a 
technique of selective electrode insulation. Namely, the 
8-in. -diameter pinch electrodes were partially insulated by 
attaching to them two circular layers of 0.005-in. mylar of 
diameter 4 in., 6 in., and 7 in., respectively (Fig. 4-9). 

In one sequence of measurements an electrode with a 6-in.- 
diameter glass center section was used to provide the desired 
electrode configuration. 

Figure 4-10 shows a series of oscillograms of driving- 
current waveforms and inner divider voltage signatures, in 
triple overlays, for the cases where 7 -in. -diameter mylar 
covered, alternately, major portions of the anode, cathode, 
and both electrodes. The quoted voltage readings were taken 
at the earliest measurable time, 0.4 jusec, when the driving 
current has a value of 130 kA and when a viable current sheet 
is becoming established at the outer wall. Figure 4-lOa shows 
the voltage signature for the normal all-metal electrode situ- 
ation. The portion of the signature of primary interest lies 
between 0.4 psec and 3.0 psec; i.e., well before the arrival 
of the current sheet at the center voltage tap. The large 
signal following at 4.0 psec is caused by magnetic flux linkage 
in the measuring circuit after the current sheet collapses at 
the chamber center. Figures 4-lOb and 4-lOc demonstrate the 
different effects on arc voltage caused by cathode and anode 
insulation, respectively. Figure 4-lOb is interpreted as 
follows: during the first 1.0 psec after breakdown the narrow 

cathode attachment has not yet contacted the electrode insu- 
lation and, therefore, the arc voltage closely approximates 
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the all-metal case over this time period. At 1.0 jasec the 
cathode attachment reaches the insulation and is arrested 
there, while the bullc of the current sheet runs inward, bend- 
ing back along the insulation to maintain contact with the 
exposed cathode. The large voltage rise between 1 psec and 
3 psec may thus be due to the increasing arc length as the 
current pattern moves inward, or to some other dissipative 
effect associated with this distortion of the cathode attach- 
ment pattern. 

When the experiment is repeated with the anode rather 
than the cathode insulated, (Fig, 4-lOc ) , the voltage is sub- 
stantially higher, 130 V compared to 70 V for cathode insula- 
tion, at the earliest observable time, 0.4 psec. It then be- 
gins to increase rapidly, indicating that the broader anode 
attachment is already being distorted by the insulating layer, 
and is thereafter arrested causing the arc length to increase 
until the sheet collapses at the chamber center. The failure 
of the arc voltage to lower to the "all-metal" voltage, even 
at the earliest observable time, is indicative of the involve- 
ment of the interior portions of the anode surface in the 
early phases of the normal pinch discharge. Once the anode 
is insulated, little change in the inception or earliest ob- 
servable arc voltage is affected by .additional insulation of 
the cathode (Pig. 4-lOd) . This result, along with that for 
cathode insulation alone (Fig. 4-lOb) , implies that the in- 
terior or central portions of the cathode do not participate, 
to the same extent as those of the anode, in the early phases 
of discharge current conduction. 

A more complete picture of these results is displayed 
in Fig. 4-11, which shows the values of the inception arc volt- 
ages, measured at 0.4 psec after breakdown, for various ex- 
posed areas of cathode and anode. The inception arc voltages 
are seen to rise from 65 +. 5 V to 135 + 5 V when the exposed 
area of the anode is decreased from 50 sq. in. (no insulation) 
to 12 sq. in. (7-in. -diameter insulation) . An equivalent de- 
crease in the cathode area causes the arc voltage to increase 
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only slightly to 80 +. 10 above data, the 

reductions in electrode area were affected by insulating the 
central portions of the anode or cathode, respectively. 

The following qualitative physical model is proposed 
to explain the anomalous arc voltage data presented in this 
section. Insulation of substantial central portions of the 
anode should not interfere with the initial formation of the 
current sheet at large chamber radii; nor should it obstruct, 
at 0.4 p.sec anyway, the development of the broad current sheet 
"anode foot" which does not become evident until later times 
in the discharge history. However, the arc voltage data of 
Fig. 4-11, corresponding to a time just 0.4 ;isec after break- 
down, does show a noticeable change in arc voltage for in- 
creasing central anode insulation. One must conclude, there- 
fore, that the central portions of the anode are also in- 
volved in the initial discharge processes. This involvement 
probably takes the form of precursor anode current conduction 
downstream of the current sheet. Although the specified anode 
insulation does not influence the initial establishment of the 
current sheet, it clearly precludes the normal development of 
precursor conduction at the central portions of the anode. 

If anode precursor conduction is an essential part of the 
pinch discharge, any obstruction of it will increase the arc 
impedance, and since the current is determined by the external 
circuit, this results in an increased arc voltage. On the 
other hand, it has been shown that similar insulation of the 
cathode has a much weaker effect on the arc voltage. There- 
fore, it is presumed that cathode precursor conduction is not 
substantial during the initial phase of the pinch event. 

These indications of anode precursor conduction in 
the large-radius pinch discharge are in qualitative agreement 
with recent es^eriments on similar high-current pulsed dis- 
charges (see Chap. II, Sec. G) . Namely, anode surface cur- 
rents, existing well ahead of the current sheet, have been de- 

29 

tected in both parallel-plate and coaxial accelerators.. In 
addition, certain unsteady plasma devices have exhibited 
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charge layers and voltage sheaths at the anode surface down- 

33 

stream of the main current distribution. Hence ^ anode pre- 
cursor phenomena, of one form or another, are common to a va- 
riety of plasma accelerators employing different electrode ge- 
ometries, current levels, and ambient gas fills. 

At this point insufficient information is available 
to warrant any attempts at proposing a mechanism for the pre- 
cursor conduction phenomena. Therefore, the collection of 
further detailed data in the central portions of the chamber 
will next be underta]<;en. 

E. ELECTROSTATIC PROBE STUDIES 

In an attempt to test the hypothesis concerning pre- 
cursor anode current conduction, various measurements were 
made of the local properties of the slightly ionized plasma 
in the region ahead of the propagating current sheet. For 
example, sensitive magnetic induction coils were inserted to 
search for possible central chamber current densities ahead 
of the main sheet. These experiments proved inconclusive be- 
cause the leakage magnetic fields due to current sheet asym- 
metries tend to overshadow those which might be related to 
possible central chamber current densities. 

Since central chamber current density measurements 
were precluded, a single tip electrostatic probe was then in- 
serted to map the plasma potential profiles in this region. 

The probe tip was spherical with a diameter of 1/32 in., and 
was formed such that it could be placed very close to the 
anode. The measuring circuit had a very high resistance, 
which forced the probe tip to float at plasma potential less 
a small sheath drop associated with the electron temperature. 

A more detailed discussion of the correction from the measured 
floating potentials to the physically more interesting plasma 
potentials is given in Chap. VI; it suffices here to state 
that the displayed floating potential distributions are a 
good representation of actual plasma potentials. 
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In its first application, the probe was used to de- 
termine the potential distribution between the pinch elec- 
trodes at a chamber radius of 1 in. with the current sheet 
standing at a radius of about 3 inches. The driving-current 
waveform used in these experiments is the same as that em- 
ployed in the arc voltage studies; the chamber was prefilled 
with 100 p argon. This distribution of potential for all- 
metal electrodes, shown in Fig. 4-12, reveals that most of 
the voltage drop occurs within 1/2 in. of the anode with a 
smaller drop at the cathode. This result relates to the 
markedly different arc voltage signatures obtained for sub- 
stantial insulation of the central portions of the anode and 
cathode, discussed in the previous section. The existence of 
a substantial voltage gradient near the anode in the region 
well ahead of the current sheet indicates that some form of 
precursor current conduction is occurring there. Hence, the 
placement of insulation on the central portion of the anode 
must inhibit this part of the current conduction pattern; the 
discharge responds to this disruption by increasing its arc 
voltage. Conversely, the near absence of precursor cathode ■ 
conduction (witness the absence of voltage gradients at the 
central cathode region) is consistent with the minimal effect 
of cathode insulation on the arc voltage. 

It is instructive to repeat the electrostatic probe 
measurements at the same 1-in. radius, but with a 7 -in. -diameter 
mylar insulation placed alternately on the anode and cathode. 
Figure 4-13 shows these potential distributions in comparison 
with that of the all-metal electrode case. Cathode insulation 
appears to have a minimal effect on the central chamber po- 
tential distribution, not an unexpected result. Anode insu- 
lation, on the other hand, gives rise to a much altered po- 
tential profile, particularly in view of the large voltage 
drop, 80 V, which appears adjacent to the insulated anode sur- 
face. This observation complies with the substantial increase 
in the arc voltage of the pinch discharge which results from 
similar insulation of the anode. The large anode voltage drop 
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is most likely caused by an electron space charge layer which 
forms over the insulator surface. The anode precursor con- 
duction, discussed on the previous page, probably involves an 
electron flux to the anode surface ahead of the current sheet. 
The presence of insulation prevents electrons from making the 
transition from the plasma to the metal anode, thus causing 
the accumulation of electron charge density there. 

The above procedures have been carried one step fur- 
ther in attempting to map the potential profiles throughout 
the 8-in. pinch chamber for the same experimental conditions. 
This was done for the case with 7 -in. -diameter mylar insulation 
placed concentrically on the anode; probe readings were again 
taken 1 psec after discharge initiation with the discharge 
current at its maximum of 310 kA and the current sheet stand- 
ing at large radii. Figure 4-14 shows the resulting eguipo- 
tential profiles. Unfortunately, their interpretation is made 
somewhat ambiguous by the fact that it is often very difficult 
to separate the true resistive components of the potential 
distribution from the inductive contributions associated with 
the unsteady features of the pinch discharge. This ambiguity, 
however, does not arise in the interpretation of the central 
chamber potential data of Figs. 4-12 and 4-13 as the inductive 
effects are minimal at these central chamber locations. Clearly, 
then, the utility of the potential profiles displayed in 
Pig’. 4-15 is considerably reduced by this problem. Neverthe- 
less, the profiles appear to indicate a concentration of po- 
tential gradients near the insulated portion of the anode with 
a much weaker field region elsewhere in the chamber. This ob- 
servation agrees qualitatively with previous results, which in- 
dicated the formation of a negative space charge layer and a 
strong potential gradient at the central region of the insu- 
lated anode. 

F. SUMMARY 

First, magnetic probe and photographic studies have 
revealed a basic asymmetry in the structure of the cylindrical 
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current sheet. Namely, the cathode region of the sheet re- 
mains narrow and well-defined, whereas near the anode surface 
the sheet becomes broad, diffuse, and sometimes divided into 
two or more regions of high-current density as it propagates 
towards the chamber center. Second, experiments using an 
inner voltage divider have demonstrated that the discharge 
arc voltage increases substantially for insulation of the 
central anode region, but that similar cathode insulation has 
only a minimal effect on arc voltage. This result has been 
interpreted .to mean that major portions of the anode surface 
particpate in the initial phase of the pinch discharge; this 
participation is presumed to be in the form of some sort of 
anode precursor conduction. The data suggests that no equiv- 
alent cathode precursor conduction exists during this initial 
phase of the discharge. These experiments were next extended 
to the parallel-plate accelerator (see Chap. V) in an attempt 
to provide a better distinction between these two anode pro- 
cesses: the spread of the current sheet at the anode pro- 

ducing the broad "anode foot" and the anode current sheet bi- 
furcation, and the apparent anode precursor conduction during 
the early phases of the discharge. 
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CHAPTER V 

EXPERIMENTS ON THE PARALLEL-PLATE ACCELERATOR 
A. PHOTOGRAPHIC STUDIES 

Studies of anode phenomena in a parallel-plate accel- 
erator have been carried out in order to determine, which of 
the effects, observed in the 8-in. pinch event, persist in 
this geometry and which are modified by it. The electrode 
dimensions in this device are 6 in. by 48 in. , substantially 
larger than those of the 8-in. pinch device. The discharge 
is driven into a 100 ja argon ambient fill by a 120 hA square 

current pulse lasting for 20 usee. In the first experiments 

41 

performed on this accelerator, a Kerr-cell photographic 
system was employed to gather qualitative data on current 
sheet structure. Figure 5-1 shows a series of such photographs 
taken through the Plexiglas sidewall of the accelerator for the 
case of all-metal electrodes. The first photograph shows an 
essentially planar current sheet standing near the back insu- 
lator wall 2 jusec after discharge initiation. After 3 psec, 
the sheet has progressed downstream and begins, to develop struc 
ture in the streamwise direction; namely, the anode attachment 
is now leading and is more broad and diffuse than the cathode 
attachment. The third picture, taken at 6 p sec, when the cur- 
rent sheet has progressed about 5 in. downstream from the back 
wall, displays the familiar "Y"-shaped structure discussed in 
Chap. IV, Sec. B. However^ in addition to the "Y" -shape of 
the current sheet, the beginnings of a "cathode foot" are also 
seen in this photograph. The final photograph, taken at 9 psec 
confirms another aspect of the anode bifurcation model proposed 
earlier; i.e., the so-called "triple point" of the Y-shaped 
sheet has now progressed towards the cathode such that the 
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anode bifurcation consumes almost the entire current sheet. 

The cathode attachment has also grown more broad and diffuse, 
a significant departure from the behavior of the 8-in. pinch 
discharge, where sharply defined cathode attachments persist 
throughout the short pinch event. This development of a sub- 
stantial cathode foot, after an elapsed time greater than 
that available to the pinch discharge, is a significant ob- 
servation to which reference will be made in Sec. C. 

B. ARC VOLTAGE MEASUREMENTS 

Inner divider voltage probe measurements, similar to 
those carried out on the 8-in. pinch device, were undertaken 
in this new geometry to further examine the importance of so- 
called precursor conduction in pulsed plasma accelerators. 

In addition, magnetic probes were inserted at various locations 
in the accelerator to permit more precise differentiation be- 
tween this precursor conduction and the broad, diffusing elec- 
trode attachments of the main current sheet. A voltage tap 
was located at the downstrecim end of the parallel-plate de- 
vice, at which position it measures only the resistive volt- 
age drop across the electrodes, exclusive of any inductive 
components. The signals were fed to an oscilloscope via a 

p 

Tektronix Model 6013 voltage probe of 10 ohm impedance. 

Figure 5-2a displays a schematic of the accelerator and volt- 
age measuring circuit; also shown is the configuration of 
electrode insulation employed in this experiment. Figure 5-2b 
shows the driving-current waveform along with a typical re- 
sponse of the inner divider voltage probe for a current sheet 
propagating along the full length of the electrodes. The arc 
voltage for this all-metal electrode case is virtually con- 
stant at a value of 80 V. 

The same measurements have been repeated with varying 
portions of the anode and cathode insulated by mylar. Typical 
triple overlayed oscillograms, corresponding to anode insula- 
tion, are shown in Fig. 5-3. Figure 5-3a, for 5 1/4 in. of 
anode exposed, shows an elevated voltage, 180 V in this case, 
immediately after breakdown followed by an increasing signal 


72 


TO 

CAPACITOR 

BANK 


W//////////ZZZA ‘ ' 7L 



CHANNEL 


’L 


MYLAR 
INSULATION 


SWITCH 



V///////////7//77/7/y7/77X V/////////7A \////7//X 


PUMP 


a) SCHEMATIC OF PARALLEL PLATE ACCELERATOR 



b) ARC VOLTAGE FOR ALL-METAL ELECTRODES 


g 

N 

(n 

m 

0. 

S PARALLEL PLATE ACCELERATOR 

CL 


FIGURE 5-2 


FIGURE 5-3 


A-3310 


A- 3315 


180- 


a) 5 1/4" ANODE LENGTH 



(HCr/»L 



— H l"*~ 2 fjL sec 


b) 10" ANODE LENGTH 


— H — 5 sec 


A- 3338 


110 - 


c) 14 1/2" ANODE LENGTH 



B 


A* 3324 



INNER DIVIDER RECORDS FOR ANODE INSULATION 






74 


as the current sheet encounters the insulation and the arc 
length begins to increase. Figures 5-3b,c,d, for increasing 
lengths of uninsulated anode, display the same general char- 
acteristics but with successively lower initial voltages and 
later onset of voltage increase. The initial elevated volt- 
ages must be again attributed to the interruption of anode 
precursor conduction by the mylar insulation. 

In order to correlate the arrival time of the current 
sheet anode attachment at the insulation interface with the 
beginning of the voltage signal increase, a magnetic induction 
coil was inserted adjacent to the anode at the location where 
insulation begins. The responses of this magnetic probe are 
shown in the upper traces of Figs. 5-3c,d for 14 1/2 in. and 
21 in. of exposed anode. Observe that the voltage signal be- 
gins to increase when the leading edge of the "anode foot" 
contacts the insulation, and then continues to rise while the 
lagging parts of the "foot" arrive there. It is not clear 
whether the voltage increase is caused by constriction of the 
previously broad and diffuse anode attachment at the insulation 
edge, or whether it is related to an increase in arc length 
resulting from the arrest of the anode attachment while the 
main portion of the sheet continues to propagate downstream. 

In either case, this voltage increase appears to be independent 
of the precursor anode conduction mechanisms which were re- 
sponsible for the elevated voltages observed at the beginning 
of the discharge. 

However, a significant departure from the behavior 
observed in the pinch event is found when the same measurements 
are repeated on the cathode. With the machine polarity re- 
versed, so that the cathode now becomes the insulated electrode, 
the arc voltage traces obtained for various lengths of cathode 
insulation (Fig. 5-4) are, in most respects, unexpectedly 
similar to those for corresponding anode insulation. Once 
again, elevated voltage readings are observed at the beginning 
of the discharge, followed by increasing voltage when the cath- 
ode attachment meets the insulation. In most cases the initial 
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arc voltage is essentially the same as that for comparable 
anode insulation, except for a moderate difference at the 
5 1/4-in. electrode situation. Initial arc voltage readings 
above the all-metal value of 80 V can only be interpreted to 
mean, for this experimental situation at least, that cathode 
precursor conduction has assumed an importance comparable to 
that of anode precursor conduction. 

The responses of the magnetic probe, located on the 

cathode surface at the insulation interface, are displayed in 

Figs. 5-4c,d along with the accompanying voltage traces. These 

signatures, in comparison with those of Figs. 5-3c,d, show 

that the current sheet cathode attachment trails the leading 

part of the "anode foot." This result confirms the so-called 

30 31 

current sheet tilt which has been reported elsewhere, ' and 
which was also observed in current sheet photographs in the 
parallel-plate accelerator (Fig. 5-1) . In addition, the mag- 
netic probe signals reveal that the cathode attachment, at 
these downstream locations, is now broad, diffuse, and even 
divided, thus displaying "cathode bifurcation. " This con- 
clusion is partially substantiated by Kerr-cell luminosity 
photographs (Pig. 5-1) . Again, the times of cathode foot 
arrival at the insulation interface correlate well with the 
first indications of a voltage increase. The arc voltage then 
continues to rise while the lagging portions of the broad 
cathode attachment reach the insulation. The cause of this 
voltage increase must be attributed, analogously to the anode 
case, to either the constriction of the broad "cathode foot" 
at the insulation edge, or to the increasing arc length caused 
by the arrest of the cathode attachment. On the other hand, 
it must be assumed that the elevated initial arc voltages are 
caused by the interruption of precursor cathode conduction 
by the mylar insulation. 

Figure 5-5 summarizes the dependence of the accelerator 
arc voltage, at 2 fisec after breakdown, on exposed electrode 
area for the two different types of electrode insulation. The 
increasing initial arc voltage for decreasing anode area is an 
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expected result, which suggests that the arc voltage again 
rises because progressively more of the anode precursor con- 
duction is obstructed by increasing amounts of anode insula- 
tion. However, the good agreement between the anode insula- 
tion curve and the cathode insulation curve represents a 
marked departure from the arc voltage data for the pinch dis- 
charge (see Fig. 4-11) . Unlike the previous situation, cathode 
insulation now has a noticeable influence on initial arc volt- 
age, thus leading to the conclusion that cathode precursor con- 
duction, for this experimental situation, has become of equal 
importance with anode precursor conduction. This discrepancy 
between the arc voltage data of the pinch and parallel-plate 
devices is taken up in the' next section. 

C. COMPAiaSON OF PARALIiEL-PLATE AND PINCH DATA 

Both photographic observations and arc voltage measure- 
ments have revealed significant differences between the various 
electrode processes occurring in the parallel-plate and 8-in. 
pinch devices. For example. Fig. 5-6, which compares the de- 
pendence of arc voltage on electrode area for these two accel- 
erators, shows that the arc voltage of the parallel-plate de- 
vice is almost equally sensitive to variations of anode or cath- 
ode area; whereas the pinch arc voltage responds much more 
severely to changes in anode area than to changes in cathode 
area. These data suggest that precursor anode conduction occurs 
in both devices, but that similar precursor current conduction 
at the cathode surface is peculiar to the parallel -plate accel- 
erator. A second, related divergence in discharge behavior is 
that a so-called "cathode foot" is evident in current sheet lu- 
minosity photographs of the parallel-plate discharge (Fig. 5-1) , 
but not in those of the pinch event (Fig. 4-5) . A well-developed 
anode foot is clearly exhibited in both cases. 

An explanation of this apparent discrepancy is sought 
by examining the- various discharge properties for these two 
experimental facilities. • First, both devices were prefilled 
to the same argon pressure and the total discharge currents were 
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approximately equal at the respective times of data collection. 
Second, although the total electrode area is somewhat greater 
in the parallel-plate accelerator, parallel-plate and pinch 
data corresponding to the same range of electrode area con- 
tinue to display a similar dependence on anode area but a no- 
ticeably different dependence on cathode area (Fig. 5-6) . How- 
ever, it should be pointed out that the pinch arc voltage mea- 
surements were tahen 0.4 jusec after discharge initiation, be- 
fore the cylindrical current sheet contacts the circular elec- 
trode insulation; parallel-plate data, not available at this 
early time because of signal noise (see Fig. 5-4) , were taken 
at 2.0 p.sec, still sufficiently early for this measurement due 
to the larger dimensions of the parallel-plate device. Thus, 
experimental indications of cathode precursor conduction in the 
parallel-plate machine, and not in the pinch device, may simply 
be a result of the difference in the times of data collection. 
This explanation would follow if the cathode precursor conduc- 
tion processes develop at a slow enough rate such that they 
are not yet observable at 0.4 jisec (pinch data), but are well 
established by 2.0 psec (parallel-plate data). Anode precursor 
conduction, on the other hand, since it is indicated by the 
0.4 psec pinch data, must develop at a faster rate. 

Further supporting evidence for this claim is obtained 
by a comparison of current sheet luminosity photographs in the 
pinch and parallel-plate devices (Figs. 4-5 and 5-1) . The 
photographs, taken at 2 jasec after discharge initiation, dis- 
play little evidence of a cathode foot in either device, although 
an anode foot is already clearly visible for both cases. Only 
at later times, say 6 jusec, does the cathode disturbance become 
well developed in the parallel-plate accelerator. Since this 
longer time period is not accessible to the pinch event, a cath- 
ode foot is never observed there. Thus, the available data 
seem to suggest that both cathode precursor conduction and the 
current sheet cathode disturbance develop at a slower rate 
than the corresponding processes at the anode. 

Lacking more detailed information at the electrode 
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regions/ the following qualitative remarks are given in ex- 
planation of the apparent time scale difference between anode 
and cathode processes in unsteady plasma accelerators. The 
regions immediately adjacent to the electrodes of a high-cur- 
rent arc discharge are characterized, respectively, by negative 
and positive space-charge zones, referred to as the anode and 
cathode falls (see Chap. II, Sec. B) . It shall be assumed, 
for the present situation, that a space-charge layer of elec- 
trons forms at the anode-current sheet interface, while a sim- 
ilar ionic space-charge layer exists at the cathode-current 
sheet interface. It is suggested that the electrons and ions 
in the respective electrode space-charge layers tend, by some 
unspecified mechanism, to spread over the electrode surfaces 
ahead of the propagating current sheet, thereby carrying the 

anode and cathode precursor conduction currents. Furthermore, 

29 

following the reasoning of Lie et al, discussed in C3iap. II, 
Sec. G, these proposed streamwise surface currents give rise 
to j X B forces near the current sheet -electrode interface, 
which are directed both forward and away from the electrode 
surface, and therefore cause the sheet to deform there. Thus, 
the development of the anode and cathode current sheet deforma- 
tions may also be related, in this way, to the proposed motion 
of electrons and ions in the downstream direction along the 
electrode surfaces. The apparent difference in the relative 
growth rates of these anode and cathode phenomena may then be 
associated with some difference between the mobilities, thermal 
velocities, or other properties of the electrons and ions in 
the respective space-charge layers. Therefore, precursor anode 
conduction and the broad current sheet anode attachments, both 
associated with fast moving electrons, should become established 
earlier in the discharge time history than similar cathode phe- 
nomena which are controlled by .much slower ions. 

D. SUMMARY 

Experiments on the parallel-plate accelerator have again 
demonstrated the existence of precursor anode conduction which 
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is distinct from, but probably not unrelated to, the broad 
and diffuse current sheet anode foot. Moreover, during the 
loi>ger time scale of this experiment both cathode precursor 
conduction and a broad cathode foot were also seen to develop. 
These latter observations led to a hypothesis wherein faster 
growth rates are assigned to the electron-controlled anode 
processes and slower growth rates to the ion-controlled cath- 
ode processes. 
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CHAPTER VI 

EXPERIMENTS ON THE QUASI-STEADY MPD ARC 
A. BACKGROUND 

The previous two chapters described experiments on 
the pinch and parallel-plate devices and pointed out several 
interesting aspects of anode phenomena, while making no at- 
tempts to explain these curious physical processes. Clearly, 
a better understanding of anode current sheet bifurcation and 
anode precursor conduction would require detailed measurement 
of the local plasma properties in the electrode regions of 
these pulsed accelerators. Unfortunately, such probing of 
the anode region plasma by electrical means was precluded by 
the severe signal noise and irreproducibility which plagued 
the data. In addition, the unsteady feature of the pinch and 
parallel-plate discharges makes data collection that much 
more difficult- In view of these obstacles, it was decided to 
conclude the experiments on the pulsed closed chamber devices, 
and to continue the study of anode phenomena in the quasi- 
steady MPD facility (see Chap. Ill) . This device has the ad- 
vantage of affording easy access to the anode region for di- 
agnostic purposes, and of providing a steady discharge phase, 
which facilitates the measurement of important plasma prop- 
erties. Furthermore, the quasi-steady MPD arc appears to be 
the most promising of contemporary plasma accelerators, and 
thus is a very appropriate facility in which to carry on these 
invest igat ions . 

First, measurements of floating potential throughout 
the arc chamber and plume of this accelerator were undertaken 
in order to determine the dominant plasma acceleration mecha- 
nisms, and to investigate various aspects of anode behavior. 
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The following pages display the resulting equipotential maps 
for a variety of operating conditions. With reference to the 
previous experiments on the pulsed devices, the influence of 
anode size on the quasi-steady MPD discharge is considered; the 
distribution of anode current density is examined in Chap. VII. 


B. FLOATING POTENTIAL MEASUEEMENTS 

The floating potential maps were constructed by means 

I 

of point -by-point probing of the discharge with a variety of 
electrostatic sensors. The probes had three configurations of 
exposed surface: (1) 1/16-in. -diameter x 1/16-in. -high cone, 

(2) 1/16-in. -diameter flat disc, (3) 1/16-in. -diameter hemi- 
sphere. It can be shown that probe size and shape are unim- 

47 

portant for the measurement of floating potential; this was 
also demonstrated empirically by comparing the readings of dif- 
ferently shaped probes, placed at the same location in the arc. 
In most cases the probes were inserted axially upstream into 
the discharge; in some cases they were inserted radially. One 
probe was bent to allow access to regions behind the anode lip. 
All leads to the probe sensing surfaces were coaxially shielded, 
and the signals were presented to the oscilloscope via a Tek- 

g 

tronix Model 6013A voltage probe of 10 ohm impedance. 

In interpreting the floating potential profiles, it is 
important to first demonstrate that they are good representa- 
tions of the more physically significant plasma potential pro- 
files. An electrostatic probe at the floating potential, by 
definition, draws zero net current, and therefore the electron 
and ion fluxes to the probe tip must be, equal. Hence, in order 
to inhibit the fast moving electrons and thus balance the elec- 
tron and ion fluxes to the probe tip, a negative Debye sheath 
forms over the probe surface. In the case of -a quiescent 

plasma, the following expression relates the floating and 

48 

plasma potent xals via the negative sheath drop 



e L ^ 177/ /J 


( 6 - 1 ) 
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where - is the nondimens ional ion current which relates 

the ion current to the probe at a given potential to the 
ion current at plasma potential. Recent spectroscopic mea- 
surements have indicated ion temperatures in excess of 10 eV 

49 

in the arc chamber of the quasi-steady MPD arc, while Lang- 
muir probe studies in this region confirm that the electron 

50 

temperatures are relatively uniform at 1.5 eV. Eurthermore, 

is taken to be about 1.8 for the given probe dimensions 
^ 48 

and plasma properties. Hence the floating potential cor- 
rection is estimated to be 


Vf - Vp « J?Te/e 


Using the value 1.5 eV for the electron temperature 
at various arc locations and operating conditions, it is 
seen that the plasma potential values are displaced some 
6 V toward anode potential from the floating values. Ion 
streaming velocities are substantial in most parts of the 
arc chamber and downstream plume, and these tend to enhance 
the ion current to the probe tip. The net effect of this 
is to decrease the magnitude of the electron retarding 
sheath, and thereby move the floating potential values closer 
to the true plasma potentials. Thus, the 6 V correction is 
conservative. In other words, the floating potential con- 
tours displayed on the following pages are in fact close 
representations of the true plasma potentials, and hence 
shall be interpreted as such. 

C. POTENTIAL CONTOURS FOR "MATCHED" OPERATION 

Recent experiments on this accelerator have been directed 
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towards the investigation of the relationship between .exhaust 
velocity and mass flow rate for constant arc currents. The 
data show that the exhaust velocity increases steadily with 
decreasing mass flow up to a point, and then remains essen- 
tially constant for further decreases in mass flow. This 
effect can be explained by assuming that the arc, at the low 
mass flow condition, is beginning to injest mass from the 
tank region and to erode the Plexiglas insulator. Thus, in 
effect, the net mass flow to the arc is no longer decreasing 
and therefore the exhaust velocity should remain constant - 
Hence, a critical mass flow is reached below which the arc 
cannot .sustain a given current level, and it must therefore 
seek additional mass through injestion and insulator erosion. 
For an arc current of 17.5 kA. these effects begin to occur 
at an argon mass flow rate of about 5.9 g/sec? 5.9 g/sec is 
thus referred to as the critical mass flow rate corresponding 
to an arc current of 17.5 kA, and when these two conditions 
coexist the arc operation is said to be "matched. " 

Figure 6-1 displays the floating equipotential con- 
tours for quasi-steady MPD arc operation at the so-called 
"matched" conditions of 17.5 kA x 5.9 g/sec. Several signifi- 
cant features of these potential contours are immediately 
evident. First, the downstream potential differs by only 
30 V (25 V if the 1-5 eV electron temperature correction is 
applied) from anode potential, and by only 15 V from the re- 
gion of gas injection. This is an order of magnitude lower 

4 

than the observed ion kinetic energy in the exhaust stream, ‘ 
and hence voids any simple collisionless electrostatic accel- 
eration model. Second, the bulk of the potential drop occurs 
within two diaimeters of the cathode and is normal to it. Note 
that this intense field region is directed so as to decelerate 
ions attempting to flow downstream, hence cannot aid in a di- 
rect electrostatic acceleration process. However, in view of 
the high current densities prevailing in the cathode region, 
along with the intense electric fields, it is clear that a 
large fraction of the total arc power is deposited in this 
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cathode envelope; therefore it must be presumed that this 
zone is primarily responsible for the high exhaust veloc- 
ities found downstream. In some portions of the envelope 

_A —i. 

the electric field may reflect a large u x B or motional 
EMF component, but in others, such as near the inlet ori- 
fices and along the axis, this component must be small, 
leaving mainly a resistive component. This, in turn, implies 
substantial ohmic heating and subsequent electrothermal con- 
tributions to the overall acceleration process. Detailed 
investigations of both the acceleration and current conduction 

processes prevailing in this cathode envelope are the subject 

50 

of another Ph.D, thesis from this laboratory. 

Turning attention to the anode region of the arc, it 
is observed that the anode fall voltage is less than 15 V 
over the upstream lip of the anode and rises to a somewhat 
higher value at the downstream face. Subsequent measurements 
of current density distribution around the anode, using mag- 
netic induction coils, have shown that most of the current 
enters the anode near the upstream portion of the lip where 
anode fall voltages seem to be lowest. This apparent de- 
pendence of anode fall on local current density has important 
implications for anode heat losses, and hence is the subject 
of a more detailed investigation (see Sec. F) . 

D. POTENTIAL CONTOURS FOR LARGE ANODE CONFIGURATION 

Earlier experiments with unsteady operation of the 
plasma pinch and parallel-plate devices (Chap, IV and V) re- 
vealed that variations of anode size had a significant effect 
on the total arc voltage. Briefly, it was found that the arc 
voltage increased substantially when the effective anode size 
was reduced so a& to obstruct the precursor anode current con- 
duction, presumed to exist downstream of the propagating cur- 
rent sheet. In view of these observations, a large 34-in. -di- 
ameter anode was fitted to the quasi-steady MPD arc, and the 
resulting arc voltage and equipotential profiles were in- 
vestigated in order to study the effects of anode size on 
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the operation of this device. Figure 6-2 displays a map 
of the floating potential contours for the familiar matched 
operation conditions, 17.5 kA x 5.9 g/sec argon, and with 
the arc chamber modified by a 34-in. -diameter anode. Note 
that the arc voltage is 160 V in this case, as compared to 
155 V for the previously discussed operation with the small 
7 1/2-in. -diameter anode (Fig. 6-1) . However, this 5 V dis- 
crepancy is most likely caused by operational deterioration 
of the tungsten cathode, and in any event is too small to 
merit further discussion. ’Furthermore, comparison of the 
overall potential distributions for similar arc operation 
with the small and large anodes reveals few noticeable dif- 
ferences, and therefore it is apparent that the availability 
of more anode surface does not influence the behavior of the 
quasi-steady MPD arc. 

Two possible explanations for this apparent disagree- 
ment with the previous experiments on the pulsed accelerators 
are immediately obvious. First, since the effects of anode 
size variation have only been observed in pulsed or unsteady 
devices, it seems reasonable to suggest that anode precursor 
conduction subsides after the transient phase of the dis- 
charge has passed, and thereafter ceases to be an important 
aspect of arc behavior, hence is not found during quasi-steady 
operation. Second, the quasi-steady MPD arc operates at 
17.5 kA, or at almost an order of magnitude lower current 
than those employed in the pinch and parallel-plate experi- 
ments, suggesting that the effects of precursor anode con- 
duction may be insignificant at the low arc current or, more 
likely, at the low ratio of current to anode area employed 
in this device. In the parallel-plate insulation experiment 

(Chap. V), this ratio, J/A, was increased to a value of 
2 

15 kA/in. before any increase in arc voltage was observed, 
whereas J/A is only 3.7 kA/in. for the quasi-steady MPD arc 
employing the small 7 1/2- in. -diameter anode. A careful ex- 
amination of the arc voltage signals corresponding to arc 
operation with the small and large anodes reveals that these 
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traces are almost identical, even during the early transient 
phase of this discharge. This suggests that the second of 
the above interpretations is the more reasonable in terms of 
present data. In any case, it is indeed clear that no basic 
improvement in quasi-steady MPD arc performance has been 
achieved by the deployment of a larger anode, and certainly 
the increase in accelerator size and weight, thereby in- 
curred, makes its use undesirable. 

Before concluding this discussion, it is noted that 
a 25 V voltage -drop (corrected to the plasma potential value) 
appears over most of the large anode surface including the 
back face (Fig. 6-2) , implying that the entire anode surface 
participates, albeit to a very small degree in some places, 
in the conduction of the total discharge current. However, 
as was pointed out previously, this small fraction of the arc 
current, which extends over most of this large anode surface, 
appears not to influence the quasi-steady operation of the 
device. 

E. VARIATION OF ARC CURRENT AND MASS FLOW RATE 

To this point, all the discussion has pertained spe- 
cifically to quasi-steady operation of the MPD facility at 
the "matched" conditions of 17.5 kA x 5.9 g/sec. Earlier 
studies have established the importance of matching the mass 
flow to the arc current for optimum operation; indeed, the 
observed characteristics were found to depend qualitatively 
on the ratio of arc current squared to mass flow, rather than 

5 

on either separately. In this section, eguipotential pro- 
files are presented for arc operation at off-matched condi- 
tions, where the ratio J /m is either greater than (underfed) 
or less than (overfed) the value corresponding to matched 
operation. 

Figure 6-3 compares the eq-uipotential contours for 
quasi-steady arc operation at three different mass flow rates, 
1.2 g/sec, 5.9 g/sec, 36.0 g/sec, and a constant current of 
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17.5 kA. In addition to a decrease in the total arc voltage 
and downstream potential as the argon mass flow rate is in- 
creased, certain differences in the potential distribution 
are readily observable. The most apparent is the develop- 
ment of large voltage gradients in the anode region for the 
1.2 g/sec mass flow, the so-called underfed case (Pig. 6-3a) . 
It is not clear at this point whether this anode envelope 
is simply a very broad anode sheath necessitated by conduc- 
tion difficulties prevailing in the lower density environ- 
ment, or whether it is in fact an increased field gradient 
of the arc column itself formed- to counteract lower plasma 
conductivities. In either event, if all of the power de- 
posited in this region were delivered to the anode the frac- 
tional anode loss at this operating condition would be quite 
high. The importance of this anode field development thus 
becomes obvious, and therefore the conduction processes par- 
ticipating in this anode envelope are investigated in 
Chap. VII. 

Turning now to the cathode region of the arc chamber, 
it is observed that a large fraction of the total arc voltage 
continues to appear across the cathode envelope, regardless 
of the mass flow rate employed, so that this region continues 
to be one of large power deposition. The importance of the 
cathode envelope for plasma acceleration can be readily dis- 
covered by noting that in both the underfed and matched cases 
the voltage drop across this region is approximately 100 V, 
compared to only 60 V in the overfed case? measurements of 
downstream ion velocities have found essentially the same 

values for the first two conditions but return a substantial 

- 4 

decrease in downstream velocity at overfed conditxons. 

Displayed in Figure 6-4 are equipotential profiles 
for three different arc currents, 42.0 kA, 17.5 kA, and 4.4 
kA, at the same argon mass flow rate of 5.9 g/sec, repre- 
senting underfed, matched, and overfed conditions respectively. 
As before the underfed condition, 42.0 kA x 5.9 g/'sec, pos- 
sesses a substantial anode envelope of voltage drop which is 
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not observed in either the matched or overfed cases. The 
voltage drop across this region is now only 25% of the total 
arc voltage; i.e., much less than the 50% fraction observed 
for the previously discussed underfed condition (17.5 kA x 1.2 
g/sec) . Since arc performance seems to depend on the ratio 
of arc current squared to mass flow rate, it follows that the 
42.0 kA X 5.9 g/sec case represents a smaller departure from 
the matched condition, and hence displays milder symptoms of 
arc underfeeding. The other basic features of the' 'quasi-steady 
arc potential distribution remain unchanged as the discharge 
.current is varied by an order of magnitude. 

F. ANODE FALL VARIATION WITH LOCAL CURRENT DENSITY 

In a previous discussion of the potential profiles 
corresponding to accelerator operation at the matched condi- 
tion, 17.5 kA X 5.9 g/sec, reference was made to the apparent 
variation of anode fall voltage with local anode current den- 
sity. Indeed the term "anode fall" has been rather loosely 
applied here because the electrostatic probes employed in 
this experiment were not designed to measure the potential 
at less than about 0.1 cm from the metal surface, and thus 
the term "anode fall" describes the total voltage drop over 
a distance 0.1 cm away from the anode surface. The discussions 
in Chap. II have pointed out that both the anode fall region 
(thickness on the order of the electron mean free path) and 
the anode adjacent £ -layer (caused by anode cooling of the 
plasma) are normally confined to a zone much smaller than 
the 0.1 cm probing distance. For the present situation, it 
is possible to approximate the thickness of the anode fall 
by computing the electron mean free path from Langmuir probe 
estimates of the electron density near the anode surface (see 
Chap. VII) . At the upstream portions of the anode the elec- 
tron mean free path is estimated to be about 0.01 cm, whereas 
at the front face of the anode this value becomes on the 
order of 1 cm. These figures probably give a rough estimate 
of the anode fall thickness at the upstream anode locations. 



96 


where ion production in the fall is important , and where 
the electron mean free path is therefore a relevant measure 
of the fall thickness. However, at the front face of the 
anode, the drift of ions away from the anode and towards the 
cathode is not significant, and thus nor is the production 
of ions in the fall. Therefore, over this portion of the 
anode surface, the fall thickness is more likely associated 
with the prevailing Debye' length, about 0,001 cm,. In any 
event, measurements of plasma potential at a distance 0.1 cm 
from the anode surface provide only an upper estimate of the 
true anode fall voltage since contributions due to possible 
£ -layer fields are included in the probe readings. How- 
ever, for convenience, the term anode fall shall henceforth 
refer to the anode voltage drop measured to within 0.1 cm 
of the metal surface. 

In order to establish the general features of the 
anode fall variation noted above, an ejqjeriment was designed 
in which the anode current density distribution was mapped 
out using magnetic probes, while anode fall voltages were si- 
multaneously monitored by electrostatic probes. These mea- 
surements were carried out during quasi-steady arc operation 
at the previously discussed matched condition, 17.5 kA x 5.9 
g/sec, as well as at two other operating conditions, 8.7 
kA X 1.2 g/sec, and 42.0 kA x 36.0 g/sec, both of which are 
considered to be matched based on their ratios of arc current 
squared to mass flow rate (see Sec. F) . Data were taken at 
three operating conditions to give more generality to the re- 
sults and to provide a broader range of anode current den- 
sities. The matched cases were chosen because the arc is 
relatively well-behaved, and the identification of the anode 
fall is less ambiguous, under these conditions. The measured 
values of floating potential adjacent to the anode surface 
were corrected by the amount ~ 4.1 kT^/e (see Eq._ 6-2), as 
it is reasonable to assume that the ion streaming energy is 
negligible over most portions of the anode; as a rough cor-- 
rection the electron temperature was taken to be constant. 
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1.5 eV, over the entire anode surface for each of the op- 
erating conditions. Figure 6-5 was then constructed by 
plotting the corrected anode fall voltages against the mea- 
sured values of local current density at the various anode 
locations and operating conditions. 

Within the accuracy of the experiments, the data 
of Fig. 6-5 appear to indicate an inverse dependence of 
anode fall on local anode current density. For purposes 
of engineering calculations it is convenient to seek a re- 
ciprocal relationship between these two quantities, since 
this would imply a constant power density at the anode fall, 

^ constant, which then could be readily inte- 
grated over the total anode area to give the power deposited 
at the anode due to the anode fall voltage. The reciprocal 
form, = 1.5 X 10 /j^ (volt/amp/cm ), plotted on Pig. 6-6, 

defines, along with the ordinate and the abscissa, an enve- 
lope within which most of the data lies; thus any calculation 
of anode power loss based on this relationship would surely 
be a pessimistic one. If the contributions due to the elec- 
tron thermal motion and the work function of the metal anode 
are included, though neglecting convection and radiation 
(see Sec. H) , a more accurate estimate of the anode power 
density for the 17.5 IcA x 5.9 g/sec operating conditions can 
be calculated, 2.0 x 10 W/cm . This figure is rela- 

tively modest in comparison to recent experimental observa- 

51 

tions. For example, Shih and Dethlefson have experimentally 

determined the anode heat flux density in argon and nitrogen 

arcs -operating at currents from 50 to 2000 A; their data 
4 6 2 

shows P^ ^ 10 to 10 W/cm ; i.e., power densities 5' to 500 

times greater than that estimated for the quasi-steady MPD 
arc. 

In conclusion, the observed inverse dependence of 
anode fall on local anode current density has major impli- 
cations for the overall efficiency of the accelerator. 

Steady flow devices operating in much lower power ranges 
are known to dissipate large portions of their input power-- 
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in some cases 60% or more — into heating of the anode. This 
dissipation is clearly related to the anode 'fall established 
by the arc, and any empirical means for reducing the fall 
voltage should be reflected in higher overall performance. 

As will be shown in Sec. H, the inverse dependence of anode 
fall on local current density results in substantially lower 
relative power loss to the anode, and hence more efficient 
accelerator operation at the higher power levels.. 

G. INTERPRETATION OP ANODE FALL DATA 

The inverse dependence of anode fall voltage on local 

anode current density, for the matched quasi-steady MPD arc 

(see Fig. 6-5) , shall now be related to the anode fall the- 

18 19 

ories of Bez and Hocker, * see Chap. II, Sec. D. Briefly, 
they distinguish between two mechanisms by which ions are 
produced in the anode fall region: field ionization and 

thermal ionization. In the field ionization case, electrons 
from the column plasma undergo essentially collisionless ac- 
celeration in the anode fall field, and then ionize neutral 
atoms upon impact. Hence, for this ion production process 
to occur the anode fall voltage must exceed the ionization 
potential of the carrier gas. Bez and Hocker have developed 
a model for this type of ionization, and thereby have de- 
rived a relation [see Eq. (2-6) ] which fits the experimental 
anode fall data for low temperature arcs, less than 7000 - 
9000 °K. For nitrogen and argon, the electron collision mean 
free path begins to decrease sharply at these temperatures 
(see Fig. 2-3), causing Eq. (2-6) to predict increased anode 
fall voltages which contradict experimental data. Bez and 
Hocker claim, therefore, that the conditions for field ioni- 
zation are violated at arc plasma temperatures above 7000 - 
9000 °K, and that ion production in the fall is then achieved 
via the thermal ionization mechanism. In this case, the ac- 
celerating electrons have their velocities randomized by 
elastic collisions in the fall zone, and since only electrons 
in the high energy "tail" of the energy distribution take 
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part in the ionization process, the anode fall voltage drops 
below the ionization potential of the gas. The important 
point of Bez and Hocker is that the transition from field to 
thermal ionization occurs at some critical arc plasma tem- 
perature, and that this critical temperature coincides with 
a sudden decrease in the electron collision mean free path 
due to increased ion-electron encounters. Furthermore, it 
was observed experimentally that the transition from P- to 

T- ionization in the carbon arc with nitrogen gas takes place 

2 

at an anode current density of 40 A/cm . 

For comparison with the above model of the anode fall 
zone, the anode fall data for the matched quasi-steady MPD 
arc have been replotted in Fig. 6-6. Also shown on this 
figure are twin Langmuir probe electron temperature measure- 
ments taken at four different anode locations for the 17.5 
kA X 5.9 g/sec operating condition. Most of the data corre- 
spending to low anode current densities, less than 60 A/cra , 
lie above the ionization potential of the working gas (argon) , 
whereas the anode fall voltages at current densities greater 
than this value are generally less than the ionization po- 
tential. In addition, the plotted electron temperature data 
indicate that the high anode fall voltages correspond to tem- 
peratures below about 19,000 °K,whereas the low voltage drops 
occur for electron temperatures above this value. Invoking 
the ionization model of the anode fall zone, it is reasonable 
to presume that anode fall voltages greater than the ioniza- 
tion potential of argon reflect the field ionization mechanism, 
while those below it correspond to the thermal ionization pro- 
cess- Consistent with this interpretation, the data define 
critical values of the near anode electron temperature and 
the anode current density: 19,000 °K and 60 A/cm^, respec- 

tively. This apparent transition anode current density is 
close to the critical value for the carbon arc, 40 A/cm^, 
thus suggesting that the quasi- steady MPD arc data may in- 
deed comply with Bez and Hocker. On the other hand, the elec- 
tron temperature of 19,000 °K is far from the argon critical 
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temperature which/ based on the sudden decrease in electron 

o 19 

collision mean free path/ is about 7500 K. Since the 
electron temperature measurements were taken at about 0.5 cm 
from the anode surface, one may argue that plasma cooling 
at the anode surface causes the plasma temperature at the 
edge of the fall zone to be substantially lower. Conceiv- 
ably, the electron temperature of 19,000 °K may correspond 
to a temperature of around 7500 °K for the cooled plasma 
adjacent to the anode fall. However, from the available 
data it is not possible to determine, one way or another, 
whether the matched quasi-steady MPD arc anode fall data 
agree, even qualitatively, with the anode fall ionization 
model of Bez and Hocker. 

H. POWER LOSSES AT THE ANODE 

The total power delivered to the anode of the quasi- 
steady MPD arc can be written as 


Pa 

- f 

cA 


!\f 5 Jkle 4- Vfi 
V a e 

da. (6-3) 




At^ODE 

Sv(KFAce 



where 

PcA 

and Pf^A 

are the 

contributions due 

to convection 


and radiation respectively, and the third term represents the 
power deposited at the anode by current-carrying electrons, 
with as the current density at the anode, Vh as the 

anode fall voltage and cf>fl as the work function of the 
anode metal (see Chap. II, Sec. E) . The term j> 5/2 (kT /e) 
represents the energy flux at the anode due to the electron 
thermal energy in the neighboring plasma? the use of the 
factor 5/2, as opposed to the factor 2 which often appears 
in this literature, is discussed in the Appendix. 
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1. Radiative Losses 

Recent experiments on the quasi-steady MPD arc have 

estimated that between 2 and 10% of the total arc input power 

52 

is lost as radiation from the highly luminous arc plasma. 

The radiative power delivered to the anode can be estimated 
by calculating the percentage of this total arc radiation 
which is absorbed by the anode surface. However, the geom- 
etry of this particular accelerator is such that only a small 
fraction of the arc radiation will be delivered to the anode, 
and by far the largest portion will be absorbed by either the 
insulating walls of the chamber or the downstream environment. 
For this reason, the term in Eq. (6-3) does not repre- 

sent a substantial anode power loss, and can be shown to be 
negligible compared with the power loss due to the electron 
current at the anode; nevertheless, the total radiative 
power loss of the accelerator may not be insignificant. 

2. Convective Losses 

Convective power losses may arise at the anode because 
of the interaction of the streaming plasma with the anode sur- 
face. An order of magnitude estimate of this contribution to 
the total anode loss is attempted by assuming that viscous 
and thermal boundary layers form at the lip region of the 
anode, and that heat is convected through the boundary layer 
to the anode surface (see Fig. 6-7) . Heat transfer by con- 
duction from the hot plasma to the remaining portions of the 
anode will be due primarily to the flux of fast moving thermal 
electrons. However, the anode current density^ ^ ^ ^ essentially 
defines this electron flux, and therefore the conductive heat 
transfer at the anode is included in Eq. (6-3) by the term 
5/2 (kT^/e) . 

The convective heat transfer at the anode is expressed 

27 

in terms of an enthalpy difference after Eckert 

Pc A - 


(6-4) 
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CONVECTIVE HEAT TRANSFER TO THE ANODE OF 
THE QUASI-STEADY MPD ARC 


FIGURE 6-7 
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where /)/ is the convective heat transfer coefficient based 
on enthalpies, is is the stagnation enthalpy of the stream- 
ing plasma at the outer edge of the boundary layer, is 

the enthalpy of the plasma in thermal equilibrium at the 
anode surface, and A,/^ is the area of the anode convection 
zone. The effects of recombination between electrons and 
ions in the boundary layer is taXen into account by writing 
the convective heat transfer in the above form. 

From downstream velocity patterns for the quasi- 

4 

steady MPD arc, it is clear that the plasma streaming veloc- 
ity at the anode lip cannot exceed 10^ m/sec. In addition, 
Langmuir probe measurements near this region of the arc indi- 
cate electron temperatures of about 20,000 and plasma 

— 9 2 

densities of approximately 7 x 10 g/cm . It is assumed 
for the present calculation that the plasma under considera- 
tion is in local thermodynamic equilibrium at a temperature 
of 20,000 °K. Argon at this temperature and at a pressure 

, ~4 53 

of 1 atm has a shear viscosity, p. 'Ph 2.2 x 10 gm/cm-sec. 
Since shear viscosity depends only -on temperature over a 
wide range of pressures and densities, this value should be 
appropriate for the present low pressure ( 0.001 atm) 

situation. Using the above estimates of plasma properties 
along with the characteristic length of the anode lip boundary 
layer, d 9^ 1. cm, the flow Reynolds number is determined 

Rey~f^sdL//C 3Z (6-5) 

This low Reynolds number suggests that the anode 
boundary layer is laminar. The convective heat transfer co- 
efficient for a laminar boundary layer on a flat plate is 

27 

given in terms of the Stanton number 

hi/fu, = Stc = . 332 / 


(6-6) 
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References 53 and 54 give the following estimates 
of the thermal conductivity and specific heat at constant 
pressure for argon at 20,000 °K and 0.1—1 atm; 

K 5.7 X 10 ^ cal/cm~sec-°K and 1-0 cal/gm-°K. 

Both of these quantities are only weakly dependent on pres- 
sure and density over a wide range, and therefore the above 
values should apply reasonably well to the low density argon 
plasma under consideration. The Prandtl number can now be 
computed 


Pr = CpA/K ~ - Ct <®-’> 

The quantities appearing in Eq. (6-6) should be 
evaluated at a reference enthalpy corresponding to the arith- 
metric mean of the plasma enthalpy at the outer edge of the 
boundary layer and the plasma enthalpy corresponding to the 
surface temperature. However, evaluation at the free stream 
conditions is sufficient for these approximate calculations. 
Substituting the above values of and . Pr , along with 

the estimated plasma density and velocity, into Eq. (6-6) 
yields the following value for the convective heat transfer 
coefficient 



0.0035 gm/cm -sec 


The stagnation enthalpy of argon plasma with a streaming ve- 
locity of lo"^ m/sec and a temperature of 20,000 °K is esti- 
mated to be, ^ 1.5 x 10 J/gm. Thus, a conservative 

estimate of the anode power loss due to convective heat trans- 
fer at the anode lip boundary, layer can now be made with the 
aid of Eq. (6-4) , 


Pcfl - hiAA(is-(A) 

» 20 kW 


(6-4) 


where the area of convective heat transfer, A , 
from the flow field diagram shown in Eig. 6-7. 


is deduced 
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The total input power to the quasi-steady MPD arc 
at the matched conditions, 17.5 kA x 5.9 g/sec, is 

= JV-r = 2.6 X lO"^ kW. Hence, the convective .power 
loss at the anode appears to be on the order of 1% of the 
input power of the accelerator, 

3. Losses Due to the Electron Current 

Since it has been demonstrated that both the radiative 
and convective heat transfer to the anode of the quasi-steady 
MPD arc are rather insignificant, the expression for the total 
power loss at the anode reduces to the form 

Pa j C'Z +Va -i- (Pa) da ( 6 _ 8 ) 

Since ~Tq and 4^^ are essentially constant over most of* the 
anode surface this ejqjression can be simplified to 

Pa ^ ^(2 ) ^- I /d ^ (6-9) 

AA/ODS. 

By restricting this formulation to the case of quasi- 
steady MPD arc operation at the three matched conditions 

sighted in Sec. F, the experimentally observed variation of 

3 2 

anode fall with current density ( ^ 1*5 x 10 W/cm ) can 

be employed here. The integral on the right-hand side of 
Eq. (6-9) can now be easily evaluated, and substituting 
T^ = 1.5 eV for average electron temperature at the anode (see 
Pig. 6-6), and = 3.5 V, the expression for the total 

anode power loss becomes 

Pfl ^ 7' 2 5 -f- 4-^ S' X 10 ^ U/ATTS (6-10) 

where J is the total arc current in amps. The fractional 
anode power loss for quasi-steady arc operation at the matched 
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conditions is then 

Pfl/Pr « 7.25/V-r + ‘^SxIO^/jVr 

where V* 7 - is the total arc voltage. 

Figure 6-8 displays this fractional anode power 
loss computed for the three matched conditions/ 42.0 kA x 36.0 
g/seC/ 17.5 kA x 5.9 g/sec, 8.7 kA x 1.2 g/sec, during quasi- 
steady arc operation. The fractional power loss drops from 
over 50% to less than 10% with an increase in arc power from 
1.5 mW to about 10.0 mW. This distribution of data dramat- 
ically reflects the observed inverse dependence of anode fall, 
on local anode current density (see Sec. F) . The implications 
of these data are clear; matched operation of the quasi-steady 
MPD arc in the megawatt range of pulse power should result in 
thermal efficiencies much higher than those attained by steady, 
low power MPD arcs which operate in the 40 to 70% range (see 
Chap. II/ Sec. E) . 
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FRACTIONAL ANODE POWER VARIATION WITH 
TOTAL ARC POWER AT MATCHED CONDITIONS 


FIGURE 6-8 
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CHAPTER Vi I 

ANODE CURRENT CONDUCTION IN A QUASI-STEADY MPD ARC 

A. ELECTRON TEMPERATURE AND DENSITY MEASUREMENTS' 

Measurements of electron temperature and electron 
number density have aided the understanding of the plasma 
potential patterns which arise in the quasi-steady MPD dis- 
charge, see Chap. VI. In addition, these electron tempera- 
ture measurements can be used to correct the previously de- 
termined floating potential values to the true plasma po- 
tentials. A description of the twin Langmuir probes and 
measurement techniques used in this experiment is given in 
detail in Chap, III, Sec. G. Briefly, these particular probes 
consisted of two separately shielded and biased 1/8-in. long 
exposed tips of 3~mil tungsten wire, oriented parallel to the 
flow axis to minimize streaming effects. This probe radius 
is less than the anticipated electron- ion mean free path in 
the regions to be probed, and larger than the Debye length. 

One of the probes was bent to allow access to the back portion 
of the anode. The probe signals were monitored with Tektronix 
P6006 voltage probes and displayed on an oscilloscope. The 
electron temperature is obtained from ratios of probe voltage 
and probe current for differently biased tips, providing both 
tips operate in the electron retarding region of the probe 
characteristic. Electron number density is determined, some- 
what less precisely, by ascertaining when the probe ceases to 
operate in the retarding region and goes over to the electron 
saturation branch of the probe characteristic. The uncer- 
tainties quoted with the following probe measurements do not 
include the effects of magnetic fields and particle collisions, 
since near the anode these should not be severe. For example. 
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the electron- ion collision mean free paths, calculated from 
temperature and density measurements near the anode lip, are 
approximately 0.01 to 0.1 cm, still at an order of magnitude 
larger than the probe radius. With the aid of magnetic probe 
measurements, the mean electron gyro radii near this part of 
the anode are estimated to be about 0.02 cm, almost a factor 
of five larger than the probe radius. 

Previous measurements with floating electrostatic 
probes revealed, for argon mass flow rates below the critical 
or "matched" level corresponding to a given arc current, that 
substantial electric fields developed in the anode region of 
the arc (see Chap. VI, Sec. E) . At higher mass flows all 
fields were confined to the cathode region of the discharge, 
except for a sheath drop immediately adjacent to the anode 
surface. To explore this effect further, a more precise and 
detailed series of experiments have been performed. Figures 
7-1, 7-2, and 7-3 display the resulting plasma potential and 
enclosed current contours for three different argon mass flow 
rates, 36 g/sec, 5.9 g/sec, and 1.2 g/sec, at a common current 
level of 17.5 kA. In addition to the previously observed 
electric field growth around the anode for the 1.2 g/sec case, 
it can be seen that a greater portion of the anode surface 
participates in the current conduction for this "underfed" 
condition. These two observations will be referred to later 
during the development of the physical model. 

Electron temperature and number density were measured 

at the locations marked on these diagrams; those near the 

cathode, labeled (1) and (2), provide reference data needed 

for the model below. Prom these measurements of electron 

temperature and number density, along with magnetic probe 

data at the same locations, characteristics of the current 

conduction mechanism in this region of the arc can be deduced. 

For example, the electron Hall parameter, , can be 

56 

estimated with the aid of the Spitzer formula for electron 
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J = 17.5 kA , rh = 5.9 g / sec 
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PLASMA POTENTIAL AND ENCLOSED 
CURRENT CONTOURS 

J = 17.5 kA » rh = 1.2 g/sec 


FIGURE 7' 3 
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collision frequency in a fully ionized gas. 


= gx 1.7/f T/yie 

(3 

- 7.Zx 10 fY)e / Te SBC 

—3 

where n is electron number density in cm and T xs the 
e e 

electron temperature in eV. Hence, 


i7e = e 


n -r- % / 

= 2.^ X 10 BTe /^£ 


(7-2) 


2 

where B is the magnetic induction field in W/m . 

Likewise, the scalar conductivity of the plasma, o~ / 

5 6 

can be computed, 

^ 53 X /t? Xn A ofin -yi 

(7-3) 

or 3 ^ 

cr = 2.7S‘ X 10 Te Mho/m 

where T is in eV. 
e 

A summary of these measurements and calculations, 
performed at the two near cathode locations labeled (1) and 
(2) for three mass flow rates, is listed in Table 7-1. The 
conduction in this region of the discharge is seen to be 
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strongly tensor in nature, based on the observations that 
SIq'^ I and ^‘/e < <r , where 
measured values of current density and electric field. This 
mode of conduction, as will be seen later, is in marked con- 
trast to the processes prevailing near the anode. 

The symbol /v in this table implies an accuracy of 
approximately 25% for the numbers quoted. 


Table 7-1 Conduction Parameters 
at Various Mass Flows (J = 17.5 kA) 


m(g/sec) 

36.0 

5.9 

1. 2 

T^{eV) 

(1) 

1.9 + 0.1 

1.6 + 0.2 

2.0 ± 

0.4 


(2) 

1.9 + 0.4 

1.4 ±0.1 

2.0 + 

0.4 

n (cm 

(1) 

14 

2.0 + 0.4x10 

1.8 + 0.6x10^"^ 

1.7 + 

14 

0.7x10 

e 

(2) 

1.4 + 0.6x10^"^ 

1.3 ± 0.4x10^^ 

1.7 ± 

14 

0.7x10-^^ 

B (w/m^) 


^ 1.4x10“^ 

^ 2.5x10“^ 


3 . 2xl0”^ 



^ 3 . 9xl0“^ 

6.1x10”^ 


5.7xl0”^ 

12 e 


0.4 

^ 0.7 


1.3 



1.8 

^ 1.8 


2.4 

0~ (mho/m) 

(1) 

^ 7.2x10^ 

5.5x10^ 


7.8x10^ 


(2) 

7 . 2x10^ 

^ 4.5x10^ 


7.8x10^ 

j/E (mho/m) 

(1) 

^ 1.1x10^ 

2.0x10^ 


1.0x10^ 


(2) 

^ 1.0x10^ 

1.5x10^ 

r\^ 

1.1x10^ 


is computed from 


B. ANODE CURRENT CONDUCTION PROCESSES 

In the anode region of the quasi-steady MPD discharge, 
both the electric and magnetic fields are generally quite weak, 
and therefore the plasma conductivity must be sufficiently high 
to enable the weak electric fields to drive the necessary 
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current. In addition, current conduction in this region may 
be augmented by diffusion-dominated electron transport to- 
wards the anode surface. Furthermore, from measurements of 
electron temperature and number density at the four marked 
anode locations, one may establish, for mass flows above the 
critical value, that the random thermal flux of electrons 
from the plasma adjacent to the anode is sufficient 'to carry 
the required anode current. Assuming a Maxwell-Bo It zmann 
distribution for the electrons, the electron current density 
at the anode -plasma boundary due to random flux is simply 

~ Ce - 0/Tle/^ AXe '] 

® ‘ 4- V 2rrfYf\e.J 

Figures 7-4a,b compare this random flux current den- 
sity, computed at the four anode locations and two different 
mass flow rates, with the actual anode current density dis- 
tribution determined from magnetic probe data. The anode 
surface has been linearized in these diagrams, so that in 
going from left to right one progresses from the insulator 
wall at the back anode around the lip to the front face. The 
favorable comparison of random electron flux with anode cur- 
rent density in these two cases verifies that the plasma in 
the vicinity of the anode is sufficiently dense and thermally 
active that the required discharge current can be passed 
through it to the anode without the need for local high elec- 
tric fields. The small potential drop right at the surface 
may be required to offset the- cooling effect of the anode 
on the plasma, and mav -'ilso serve to capture those electrons 
which approach the anode by virtue of their random motion. 

One may even speculate that it is this match of ambient plasma 
conditions to the discharge current that identifies the "crit- 
jcal" mass flow for a given arc current operation. This spe- 
cification of the critical mass flow for a given arc current 
defines the so-called "matched" operation referred to in 
Chap. VI. 



j A , A/ cm^ i , A/cm 



ANODE BACK LIP ' ANODE FACE 


( Q ) m = 36.0 g/sec 



ANODE BACK LIP ANODE FACE 


( b ) m = 5.9 g /sec 

COMPUTED ELECTRON RANDOM FLUX COMPARED TC 
MEASURED CURRENT DENSITY AT ANODE SURFACE 


FIGURE 7-4 
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Regrettably, it has not been possible to make the 
same comparison for the "starved" condition, m = 1.2 g/sec. 

In this case the electrostatic and magnetic probe data in 
the vicinity of the anode are exceedingly noisy and irre- 
producible, indicative of an unhappy domain of arc operation, 
and it is therefore impossible to define characteristic 
values of the plasma properties and current densities over 
most of the anode surface other than recalling that, at this 
condition, a sizable electric field extends far from the 
anode surface, even to the edge of the cathode region (Fig. 
7-3) . In the absence of better data at the anode, the fol- 
lowing hypothesis is advanced. At the low mass flow rate, 
the electron number density near the anode is lower than 
that demanded by the random flux condition invoked above, 
presximably due to the lower overall density level in the 
chamber. Hence, in order to satisfy the current demands at 
the anode, the arc establishes an anode region field whose 
function is to increase the electron flux at the anode to 
the proper level. In addition, at this condition the arc 
current spreads itself over a larger portion of the anode 
surface (see Fig. 7-3) , thereby alleviating somewhat the 
local demand for electron random flux. 

It is possible to distinguish between two means by 
which this anode region electric field may accomplish the 
needed increase in electron flux to the anode surface: a) 

it may heat the plasma locally by means of the Joule-heating 
component of the total E * j energy input, thereby increasing 
the ionization level and possibly also the electron tempera- 
ture to the point where the random fltox of electrons again 
becomes adequate; b) it may establish an electron conduction 
current of the required magnitude to just compensate for the 
deficiency of random electron flux at this low density situ- 
ation. Two elementary calculations should serve to estimate 
the relative effectiveness of these two means. 

First, the complete expression for the electron flux 
to the anode is derived by considering a Maxwell -Boltzmann 
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distribution that has been shifted by the drift velocity 
of electrons towards the anode, Ue $ and by integrating 
over one-half of the velocity space? i.e., over all those 
electrons moving towards the anode surface. Taking the 
x-direction to be perpendicular to the anode surface this 
has the form: 



( 7 - 5 ) 

' * [(f#)1] 

where c^ = (8 kT^/rrm^)"^ is the mean thermal velocity of the 
electrons. For the high temperature plasmas that arise in 
discharges of this nature, the electron thermal velocities 
are much greater than their drift velocities, hence 

ai c< 2 xj^"na/pntf- and u.<^ «Ce. , and the lowest ap- 
proximation to ia just the expression used earlier 

for the isotropic electron distribution. 


= mee cg e 
4 




-f- 


* 

jed « ^'r.k = errie^/'f 

Now, let and define the in- 

creases in electron flux to the anode produced by Joule-heat- 
ing of the plasma and by electron conduction respectively, 
where 



'b 4ef\ \ A/T)e 



( 7 - 7 ) 




A Ha 


( 7 - 8 ) 
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The relative importance of the two terms in the ex- 
pression for depends strongly on the state of 

the plasma under consideration. For example, if the plasma 
is only partially ionized, further energy input will increase 
the degree of ionization with the electron temperature 
charging only slightly until full ionization is approached. 
The addition of energy to a plasma which is almost completely 
ionized, on the other hand, will be reflected by an increase 
in the electron temperature rather than by the onset of sec- 
ond ionization because of the high-energy threshold of the 
latter process. Unfortunately, few data are presently avail- 
able on the ionization levels in the anode region of the 
quasi-steady MPD discharge. However, measurements of elec- 
tron temperature at the back surface of the anode for the 

three different mass flows have revealed that T does not 

e 

increase at the "starved" case in response to anode field 
development . 


Furthermore, partial differentiation of the expression 
for the electron flux to the anode, Eq. (7-5) , yields the 
following results: 


A(7)q 

'bma) 


J 


ATe 


COqQCq 



z 


/\Cr\e. 


-2. 

— U.Q. 

(Tiee^ e ^1^ 

'I- a7^ 


(7-9) 


(7-10) 


where the assumptxon, Ue « 2 ^ le has again been employed. 
These expressions show that the fractional increase in elec- 
tron temperature must be double that of electron number den- 
sity in order to achieve the same increase in the electron 
flux to the anode. Hence, it becomes apparent from both em- 
pirical and theoretical arguments that an increase in the level 
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of ionization/ rather than a raising of the electron tem- 
perature/ is the more likely and more effective means by 
which Joule-heating of 'the anode region plasma may serve to 
raise the fl\ax of electrons to the anode. Therefore the 
second term in the expression for (A may be neglected; 

and hence a comparison of the relative importance of plasma- 
heating and electron conduction involves the evaluation of 
the two partial derivatives: and 

With reference to Eq. (7-5) , the following expressions 
can be derived; 

a 

- /me 

= g ~^e e (7-11) 

^ (vx^e I / + 2 ^ 

Employing the usual assumption that Ltg « 
to simplify the above forms / a ratio is established which com- 
pares the relative effectiveness of plasma -heating and elec- 
tron conduction in terms of their ability to raise the elec- 
tron flux at the anode to the required level/ namely 

— f 2 ^ ^ /SCCiQ. _ 

(A{ef^)c me/iUe 

where is the normalized increment in electron density/ 

and AUe is the increment in electron drift velocity per- 
pendicular to the anode normalized by the electron mean ther- 
mal velocity. 

The respective increases in electron number density 
and electron drift velocity/ and AUq , which arise 


^(V)e 


(7-13) 


VytT 

“b jeA 
'6Ue 
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in response to the anode region field at the starved con- 
dition, , must now be related to this field by means of 

ionization and conduction models. First, the current con- 
duction in the anode region of the discharge is probably 
tensor in nature,- with the components of current parallel 
and perpendicular to the anode field given by ) 

and i perpendicular com- 

ponent lies in the r~z plane and is directed towards the down- 
stream portions of the anode. Hence, the increase in electron 
drift velocity perpendicular to the anode surface, i.e., par- 
allel to E*a / is given by 


A U e ^ ( -A /y X (Y ) e £* ) (5~ ^ Ea / ( / -f-_0 e } 




tnTA\ 


where is the increase in the anode region field, which 

results because of the "starved" mass flow rate, and the ratio 
<y/{r\e is assumed to be only weakly dependent on /TI0 

The rate of Joule-heating in the . anode field region 
2 57 

is simply j /cr , where j is the total current density given 
by' ^ • Thus, the rate of ionization, or 

in other words, the rate of electron production due to the in- 
creased anode field, Ept , is given by 


d(Y\e. 

dt 


/Vcr ^ q-ETa 
e V/ eVi* 


(7-15) 


where it is assumed that all of the Joule-heating energy goes 
into ionization and i® the. effective or mean ionization 

potential of argon. In order for the steady state level of 
electron density at the anode to increase by the amount 
over the value it would assume in the absence of ionization, 
^rt^^new electrons must be produced in the anode field region 
over the time interval taken for all of the electrons in that 
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region at a given moment to reach the anode. In other words, 
the time interval, , used to compute / must be 

approximately equal to the transit time for electrons drift- 
ing across the anode field region. If the average thickness 
of this region is , then the time interval during which 

electrons are produced must be = cLa /u<s. ^ 

where Ug is again related to the current component par- 
allel to £a / ///\<£ e • Hence the increase 

in electron number density at the anode, due to Joule-heating 
of the plasma in the anode field region, is given by 

^ — g » cjfK (7-16) 

ev* ere'// Ue 


which reduces to the simple form 


A/r?e 


Ale d / \/i 




(7-17) 


Substituting the expressions for and A//)e given 

by Eqs. (7-14) and (7-17) into the ratio of Joule-heating to 
electron conduction defined by Eq. (7-13) yields 

fA^efl)cr _ /2 JjTe /ng e (/+n|) 

erV* 

Using the formula for (5~ given by Eq. (7-3) , and 
substituting cLa ^ 2 Cy] for the approximately uniform thick- 
ness of the anode field region^ and = 33 V for the effective 

5 8 

or mean ionization potential of. argon, the ratio takes the 
simplified form, 

(^A hO X (0 tY\s C I + -O-e} / 

(a 

O 

where me is in cm and Te is in eV. 


(7-19) 
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The above ratio is strongly dependent upon local 
plasma properties; hence must be evaluated at several re- 
presentative anode locations in order to determine whether 
Joule-heating or electron conduction is the more effective 
means of increasing the electron flux to the anode surface. 

As mentioned earlier, considerable noise and irreproducibility 
are associated with Langmuir probe signals at this "starved" 
condition, thus precluding the measurement of electron tem- 
perature and number density in the anode region of the arc. 
Lacking better data, the value T" e =2 eV (see Table 7-1) 
shall be assumed for the plasma at the lip region of the 

anode. Measurements of electron density at this same loca- 

14 —3 

tion for the 5.9 g/sec mass flow yield fY\e. ^ 1 x 10 cm 

(see Fig. 7-4b) . Presumably, the electron densities are 
lower at the starved case because of the lower overall plasma 
density at - the decreased mass flow, and as a first approxi- 
mation, it is assumed that the electron densities are lowered 

by the ratio of the respective mass flow rates; hence 
13 — 3 

me ^ 2 X 10 cm is predicted at the anode lip for the 

1.2 g/sec mass flow case. The magnetic fields at the lip 
region are estimated from the enclosed current contours of 
Fig. 7-3, thus 6 ^ 3.5 x 10 W/m^. However, when these 

values of ~e , fAe. , and 8 are inserted into the formula 
for the electron Hall parameter, 

'5 5/p / 

-Oe = X iO B Te / 


the result is J2e ~ 12. A Hall parameter of this magnitude 
implies that the current near the anode must flow essentially 
perpendicular to the electric field. This is clearly in dis- 
agreement with the current and potential contours for the 
starved condition (Pig. 7-3) , suggesting that the first order 
approximation for nOe based on the ratio of mass flow rates 
is not valid. Magnetic probe data, used in constructing the 
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enclosed current contours for the starved case, were not 
available at less than about 1 cm from the anode surface 
because of excessive signal noise there, and hence the cur- 
rent lines close to the anode are very approximate. Never- 
theless, by comparing the direction of current flow near the 
anode surface, as indicated by the enclosed current contours, 
with the direction of the anode region field. Hall parcimeters 
between 0.5 and 2 are estimated over much of the anode region. 

Looking first at the lip region of the anode where 

14 —3 -T- 

rrig should be approximately 10 cm and is about 

2 eV, the ratio of Joule-heating to electron conduction be- 
comes 

{a ^ 

(Aje^)cJ 

/\a/00&.LIP 

Thus, it appears that I 

sition regime between scalar and Hall conduction, also im- 
plies a condition at the anode lip wherein the Joule-heating 
and electron conduction mechanisms assume aJLmost equal im- 
portance in maintaining the required electron flux to the 
anode surface at the starved condition. 

At the front face of the anode, the magnetic fields 
drop off by factors of at least two to seven, depending on 
the distance from the anode lip; since SXe is relatively 
constant over most of the region in question (see Fig. 7-3) , 
Eq. (7-2) requires that the electron density must decrease 
by similar amounts. The electron temperature may also de- 
crease at the front anode but not sufficiently to compensate 
for the substantial decrease in * therefore the ratio 

by Sg. ( 7 - 19 ), yields at 

the front face of the anode. 




(7-20) 


, which characterizes the tran- 


1 

(a ^'e A )c j 

/^vooe FACE 
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This result implies, for this low density, low magnetic field 
region of the anode, that an increase in the electron flux 
to the anode surface at mass starvation is more effectively 
accomplished by the process of electron conduction than by 
that of Joule -heating. This conclusion is also physically 
intuitive since the low plasma densities tend to discriminate 
against the ionization mechanism, while the weak magnetic 
fields encourage current conduction in a direction parallel 
to the anode field, i.e., towards the anode surface. . 

An analogy between the quasi- steady MPD arc anode, 

whose current demands are satisfied simply by random thermal 

electron flux from the adjacent arc plasma, and a Langmuir 

probe, which is operating at the electron saturation point 

of the probe characteristic, is suggested. Several elaborate 

collisionless theories have been developed to describe the 

sheath growth around a probe tip which is forced to draw 

more than the saturation, or random thermal, electron cur- 
44 

rent.' One may speculate that the anode field growth ob- 
served at the starved mass flow, or low plasma density, con- 
dition is in fact an electron capture sheath similar to that 
formed around an electron attracting probe tip. Should this 
be true, Langmuir probe theories could be used to describe 
the anode field region of the starved quasi-steady MPD arc. 
However, calculations of Debye length and collision mean free 
paths in the MPD arc plasma reveal that both are several 
orders of magnitude smaller than the anode field envelope. 
Hence, the process of anode field growth at mass starvation 
is in no way similar to the collisionless Debye sheath which 
grows around Langmuir probe tips attempting to draw more than 
the saturation electron current. 

C. SUMMARY 

It has been shown, for the quasi-steady MPD arc at 
matched and overfed mass flows, that anode current conduction 
is accomplished by the random thermal flux of electrons from 
the adjacent arc plasma. At conditions of mass starvation. 
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substantial electric fields develop around the anode. It is 
proposed that the random thermal electron flux is insufficient, 
at this low mass flow, low plasma density situation, to carry 
the necessary anode current, and hence these anode region 
fields are formed to maintain the electron flux to the anode 
surface at values consistent with the arc current demands. 

In the anode lip region of the arc, the processes of electron 
tensor conduction and Joule-heating of the plasma both con- 
tribute to this task, whereas at the front face of the anode, 
the electron conduction mechanism appears to dominate. 
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CHAPTER VIII 
SUMMARY AND CONCLUSIONS 


A study of anode phenomena in high-current arc dis- 
charges has been carried out under several different es^eri- 
mental situations employing various electrode geometries, arc 
current levels, and current pulse durations. The unifying as- 
pect of these experiments, on several rather different plasma 
accelerators, has been the -common goal of achieving more ef- 
fective plasma acceleration through the understanding of the 
various physical processes occurring at the anode of these 
devices. Specifically, the pulsed or unsteady mode of gas ac- 
celeration was studied in both cylindrical pinch and parallel- 
plate geometries at current levels from 120 to 300 kA, pulse 
durations from 3.5 to 20 psec, and ambient argon pressures of 
100 p. The quasi-steady mode of accelerator operation, in a 
coaxial MPD arc geometry at current levels from 4.4 to 42.0 kA 
and argon mass flows from 1.2 to 36.0 g/sec, was also examined 
in detail. 

First, investigations of current sheet behavior in the 
large-radius pinch discharge and parallel-plate accelerator 
have revealed that the narrow current distribution begins to 
deform at the anode surface into a broad, diffuse, and some- 
times divided "anode foot." At later times, a similar struc- 
ture, though somewhat less pronounced, develops at the cathode, 
and this "cathode foot" lags behind the anode attachment as the 
current sheet propagates inward or along the electrode channel. 
A random, filamentary "current spoke" has been '.observed at the 
anode of the pinch discharge, and combined with irreproducible 
probe signals in this part of the chamber, indicates an uncer- 
tain or unstable domain of current conduction. A similar con- 
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duct ion difficulty is encountered, at underfed mass flows, in 
the anode region of the quasi-steady MPD arc, as indicated by 
excessive noise and irreproducibility of magnetic and electric 
probe data. These observations may also imply substantial 
power dissipation in the anode region of both the pinch and 
MPD discharges. 

In addition to the main current distribution in the 
pulsed accelerators, a less obvious but equally important re- 
gion of current conduction appears to form at the electrode 
surfaces considerably downstream of the current sheet. The 
pinch discharge experiments show that this precursor conduc- 
tion occurs mainly at the anode surface, but the longer time 
scale, parallel-plate experiments reveal that the cathode also 
participates in this downstream current conduction, although 
to a lesser degree. A qualitative argument, based on the 
relative mobilities of electrons and ions, was presented to 
expiain the slower growth rates of cathode conduction processes 
compared to those at the anode. However, the real significance 
of precursor current conduction is its apparent relation to the 
power dissipation of pulsed accelerators. More explicitly, the 
resistive arc voltage of both the pinch and parallel-plate de- 
vices was found to increase sharply when precursor anode con- 
duction was impeded by mylar insulation; similar interruption 
of precursor cathode conduction in the parallel-plate device 
induced a somewhat less severe increase in resistive voltage 
drop. From results of this nature, it can be proposed that 
the thermal efficiencies of pulsed plasma accelerators should 
decrease with a decrease in electrode area, with the effect 
being more strongly dependent on a decrease in anode area. A 
very different result was discovered from experiments on the 
quasi-steady MPD arc wherein the anode fall voltage tends to 
decrease with increasing local anode current density, thus 
implying that higher thermal efficiencies may, in fact, be 
achieved by decreasing the anode size. 

Floating potential measurements on the 8-in. pinch dis- 
charge revealed that almost all of the resistive arc voltage 
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drop occurs within about 1/2 in, of the anode surface, with 
a sizable portion of this voltage appearing adjacent to the 
surface. These data imply that most of the power dissipation 
in this device takes place near the anode; this observation 
carries over to the quasi-steady MPD arc, where detailed po- 
tential patterns indicate that the anode is a major source of 
power loss, while most of the actual plasma acceleration takes 
place in the cathode region of the arc. 

As pointed out in Chap. I, high power quasi-steady 
MPD arcs appear to promise high efficiency, high specific im- 
pulse, high thrust density plasma acceleration; certain re- 
sults of- the present study speak directly to the high thermal 
efficiency of these devices. For example, an apparent inverse 
dependence of anode fall voltage on local anode current density 
was observed for matched quasi-steady operation of the MPD fa- 
cility at current levels of 8.7, 17.5, and 42.0 kA. When the 
empirical relationship, fig 1.5 x 10^/j^ (mks) , was used in 
a computation of the fractional anode power loss, the results 
predict a decrease from over 50% to less than 10% when the total 
arc power is increased. from 1 to 10 mW.. These calculations 
were based on -the assertion that the electron current is the 
major cause of anode heat transfer, whereas convective and 
radiative losses at the anode are much less significant. 

Finally, the importance' of mass flow rate to the over- 
all characteristics of the quasi-steady MPD discharge, in par- 
ticular to the arc behavior in the anode region, was determined. 
Namely, at matched and overfed argon mass flows, the electric 
fields, except for those in the acceleration region near the 
cathode, are confined to a narrow region at , the anode s^arface; 
and anode current conduction is accomplished simply by random 
thermal electron flux from the adjacent arc plasma. However, 
at underfed mass flows substantial fields develop around the 
' anode. Approximate calculations indicate that these fields 
help to maintain the required electron flux to the lip region 
of the anode by establishing both Joule-heating of the plasma 
and tensor electron conduction; at the front face of the anode. 
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electron conduction appears to be the dominant mechanism for 
producing the requisite electron flux. The power deposited 
in this anode field envelope probably appears as an increase 
in fractional anode power loss; in any event, it is respon- 
sible for a decrease in the thermal efficiency of the accel- 
erator at the underfed condition. Hence, the available data 
imply that high-efficiency plasma acceleration can be accom- 
plished by quasi-steady MPD operation in the megawatt range 
of pulse power, provided that sufficiently high propellant 
mass flow is supplied to avoid the formation of electric 
fields in the anode region of the arc. 
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Appendix 

CONTRIBUTION OF ELECTRON THERMAL ENERGY 
TO ANODE HEAT TRANSFER 

For some time there have been conflicting opinions re- 
garding the magnitude of one of the terms appearing in the 

anode heat transfer expression generally applied to MPD arc- 
5 9"^5 2 X3 

jets. ' The expression is usually written in the 

following form, 

= °RA + °CA J(5kT^/2e + (A-1) 

where is total heat transferred to the anode and are 

the heat transferred by radiation and convection. The third 
term represents the energy carried to the anode by the elec- 
trons which are presumed to carry most of the current and 
is anode fall voltage, is work function of metallic anode, 
and T^ is electron temperature. 

In question is the numerical coefficient of the first 
term in the bracket which represents the thermal energy trans- 
ported to the anode, per current-carrying electron. Some au- 
thors prefer the value 2kT^/e, arguing that the electron pres- 
sure contributes only l/2kT^ to the enthalpy term because cnlyone-. 
half of the electron distribution is absorbed by the anode. 

Hence 

^ = f f 2 <^-2) 

= - (3/2kT + l/2kT ) 

e e • e 

= J (2kT /e) 
e 
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where h is electron specific enthalpy and P , 6 are electron 
e e e 

pressure and density. 

We shall now look at this problem from a particle flux 

point of view and discover that both of these forms have some 

supporting evidence but that the stronger arguments point to 

5kT /2e as the correct form, 
e 

First, let us assume that the electrons have a normal 
Maxwell -Boltzmann distribution at the edge of the anode sheath, 
and that all electrons having velocities directed towards the 
anode are captured by the voltage drop of the fall and absorbed 
by the anode. Defining the positive x-direction as the inward 
normal to the anode surface, we can write the energy per unit 
area carried to the anode by electrons: 
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With reference to standard integral tables the result 


xs 


irn m 

T = — ^ ^ 

*Ae 2 


/2kT 

(^) 


(A-4) 


The flux of electrons to the anode is given by a similar 
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integral 
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Using this result in our previous expression we then find that 

■SAe “ (A-6) 


which can be integrated over the anode area to give 

Q^e " J(2kT^e) (A-7) 


In the second approach we employ a Maxwellian distri- 
bution for the electrons that has been shifted by the electron 
drift velocity so that there is now a net average motion of 
electrons towards the anode. The particle integration is 
carried out at a surface near the anode sheath edge but never- 
theless in a region where electrons can cross the surface in 
both directions. If one assumes that there is no electron 
heat flux through the side boundaries of the arc then the heat 
flux across the surface equals that entering the anode. The 
positive x-direction is defined in the same way as before, but 
the integration is now carried out over all particles. The 
electron heat flux is now given by 

/ v3/2 ^ «o 

“SAe = -f ”e(2^) / / / V ^ 


[ (V, 


— *^e/ e 


u) 


+ v. 


+ 


dv dv dv 
X y z 


(A-8) 



136 




o m /2kT (v - u) 
3 — e' e x 


dv + kT n u 
X e e 


After making the change of variables 13 = v — n the integral 
transforms to ^ 
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Again referring to the tables and employing symmetry 
arguments;, the final form becomes 


= 3/2kT n u + l/2n m u"^ + kT n u 
^Ae e e e e e e 


kT m u 
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By physical reasoning or by carrying out the integration 
for electron flux density across the imaginary surface 


oo 


00 " cO — oo 


V dn = n u 
X e e 


(A-11) 


Hence our expression for the electron energy flux density 
becomes. 


‘^Ae = 3^5/2 -f ^1/2 -f 


(A-12) 


Realizing that electron thermal energies are much 
greater than their drift energy for our plasma, and integrating 
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this expression over the surface area we arrive at 

°Ae = JC5kT^/2e) (A- 13) 

which is the most frequently quoted form of the electron energy 
flux to the anode. 

As derived, the difference between the two calculated 
forms of is seen to arise from the presumption of the first 

method that an isotropic Maxwellian distribution exists at the 
edge of the anode sheath. In view of the steady demand for de- 
livery of electrons of one velocity sign to the sheath, this- 
assumption is untenable in this geometry in the steady state. 

To maintain such a condition, the sheath edge would need to 
recede from the anode surface as time progressed. The second 
method makes no such presumption, hence yields the proper re- 
sult. 
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